Black Hole in the Expanding Universe from Intersecting Branes 
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We study physical properties and global structures of a time-dependent, spherically symmetric so- 
lution obtained via the dimensional reduction of intersecting M-branes. We find that the spacetime 
describes a maximally charged black hole which asymptotically tends to the Friedmann-Lemaitre- 
Robertson- Walker (FLRW) universe filled by a stiff matter. The metric solves the field equations 
of Einstein-Maxwell-dilaton system, in which four Abelian gauge fields couple to the dilation with 
different coupling constants. The spacetime satisfies the dominant energy condition and is character- 
ized by two parameters, Q and r, related to the Maxwell charge and the relative ratio of black-hole 
horizon radii, respectively. In spite of the nontrivial time-dependence of the metric, it turns out 
that the black hole event horizon is a Killing horizon. This unexpected symmetry may be ascribed 
to the fact that the 11-dimensional (11D) brane configurations are supersymmetric in the static 
limit. Finally, combining with laws of trapping horizon, we discuss the thermodynamic properties 
of the black hole. It is shown that the horizon possesses a nonvanishing temperature, contrary to 
the extremal Reissner-Nordstrom (RN) solution. 
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I. INTRODUCTION 

Brane configurations in string/M-theory have 
offered a new avenue for producing wide classes 
of solutions of physical interest in lower dimen- 
sions. The attempts of an early date have mainly 
aimed at the construction of various kinds of black 
holes [ll-Q via the Kaluza-Klein compactification 
of (intersecting) branes An interesting ap- 

plication of this idea is to dynamically realize our 
4D universe by incorporating time-dependence fiol — 
Il2| . An alternate mechanism that provides lower- 
dimensional spacetimes is the brane world [l3|, [l4| , 
in which our 4D world is regarded as the hyper- 
surface embedded in the bulk. Cosmological evolu- 
tionary scenarios based on the dynamically moving 
brane or colliding branes give a significant modifi- 
cation from the standard cosmology in high energy 
regime, still consistent with the present day obser- 
vations [TB-fiS. 



One can extend these studies further into the 
case where the brane involves a nontrivial space- 
and time-coordinate dependence. Correspond- 
ingly, the 4D reduced solution becomes spatially 
inhomogeneous and evolving in time. A prelimi- 
nary result was presented in [l8| . where colliding 
D3 branes were discussed within the framework of 
type IIB supergravity (see [l9[ for analysis of the 
Hofawa-Witten domain wall). Lower-dimensional 
solutions obtained by compactifying extended di- 
rections of these moving branes or dynamically in- 
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tersecting branes have much richer properties than 
those obtained from static counterparts. Recently, 
the authors of [2(| have studied dynamical solu- 
tions describing intersecting branes in more gen- 
eral settings and obtained a number of intriguing 
solutions with wide potential applications. Among 
other things, their tantalizing findings are the "cos- 
mological black holes" in the expanding FLRW uni- 
verse, that is, black holes in a non-isolated sys- 
tem which asymptotically approaches to the ho- 
mogeneous and isotropic cosmology. In this paper, 
we are concerned with this "candidate" black hole 
spacetime obtained in f20j . 

Studies of black holes in our universe have been 
primarily focused upon the stationary spacetimes 
in the literature |2l| . Such approaches are defi- 
nitely the first step because we can anticipate that 
dynamic variations will die away and the system 
will settle down to equilibrium states if sufficiently 
long time passed after the formation of a black hole. 
A number of physical properties of stationary black 
holes have been elucidated by many people. Specif- 
ically, black hole uniqueness theorem in stationary 
spacetime is the major triumph of mathematical 
relativity and established that the Kerr solution 
describes all vacuum black holes in isolated sys- 
tems [2l| . An essential crux toward the uniqueness 
proof is the demonstration that the event horizon 
in stationary spacetime is a Killing horizon (22j . 
Since a Killing horizon is a null surface to which 
the Killing vector is normal, we can identify lo- 
cus of a black hole simply by the local spacetime 
symmetry. Furthermore, it has been revealed that 
Killing horizons admit three laws of black hole me- 
chanics which bear an amazing resemblance to or- 
dinary thermodynamics [23l430| . This implies a 
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deep association between classical gravity, statisti- 
cal physics and quantum mechanics, so that black 
hole thermodynamics is expected to have a key role 
for understanding quantum aspects of gravity. A 
notable progression in string theory is the micro- 
scopic derivation of black hole entropy in the 
spectivc of intersecting brane configurations 

If we get rid of the stationarity assumptions to 
discuss dynamics, uniqueness theorem no longer 
holds. Accordingly, a variety of black hole solu- 
tions are likely to exist. However, very little has 
been known concerning the exact solutions of Ein- 
stein's equations that describe growing black holes 
interplaying with surroundings. A novel aspect of 
non-isolated system is that it generally possesses a 
time-dependence and need not be asymptotically 
flat. The background fluid distributions filling the 
universe become inhomogencous due to the pres- 
ence of the black hole, while the black hole grows 
by swallowing the ambient matters. Such a com- 
plexity has rendered the system considerably elu- 
sive. 

A large amount of effort has been devoted thus 
far to attempt to obtain black holes in the FLRW 
universe. An initiated work is the simplest model 
invented by Einstein and Straus |3l[, which is of- 
ten refereed to as a "Swiss-Cheese Universe." They 
matched the Schwarzschild solution with an FLRW 
universe by means of a "cut and paste" method. So 
the black-hole part metric still maintains a time- 
symmetry, and then it appears that this model 
does not capture realistic situation of dynamic 
phase. If there is a positive cosmological con- 
stant, we have the Schwarzschild-dc Sitter (SdS) 
or the Reissner-Nordstrom-de Sitter (RNdS) solu- 
tion [U HH . Those spacetimc can be redescribed 
by the coordinate transformation in the form of 
a black hole in the exponentially expanding uni- 
verse [H|]. However they have a "timelike" Killing 
vector and are essentially static. 

In recent years, Sultana and Dyer have con- 
structed a more sophisticated black hole solution 
in a dynamical background by a conformal tech- 
nique 3J] . The matter content is composed of null 



dust and usual dust fluids, and the solution tends 
to an Einstcin-dc Sitter spacetime asymptotically. 
This model, however, suffers from the issue of vio- 
lating energy conditions: energy densities of both 
fluids become negative at late times. 

One of the other widely known black hole can- 
didate in FLRW universe is the McVittie solu- 
tion [35j . which is a spherically symmetric, expand- 
ing solution of the Einstein equations sourced by a 
perfect fluid. Taking asymptotic limits, the solu- 
tion looks like an FLRW universe at "infinity," and 
like a black hole near the "horizon," whence one 
might be led to conclude that the McVittie space- 
time might describe a black hole in the expanding 
FLRW universe. Attractive as this might be, how- 



ever, such an optimistic outlook would jump to a 
hasty conclusion. As asserted in (36|, the McVit- 
tie solution is disqualified as a black hole in an 
FLRW universe. Since our spacetimc metric is in 
appearance quite similar to the McVittie solution 
in several respects, the above concrete example mo- 
tivates us to explore the global structure of time- 
dependent black holes with enough care. In this 
paper, we intend to provide a comprehensive ac- 
count of the global picture of dynamical solutions 
obtained in [2G |. 

Another interesting issue of a time-dependent 
spacetime containing black holes is the collision of 
black holes. Kastor and Traschcn found the col- 
lection of extremely charged black holes in the de 
Sitter universe and discussed their collision [37], H|| . 
This solution is a time-dependent generalization of 
its celebrated cousin, the Majumbdar-Papapetrou 
solution, which describes extremely charged RN 
black holes . In spite of the lack of Bogomol'nyi- 
Prasad-Sommerficld (BPS) states for A > 0, it 
is somewhat astonishing that the superposition of 
RNdS black holes is possible. The same situation 
happens in our spacetimc obtained by the time- 
dependent intersecting brane system in which no 
BPS states are preserved. Making use of this exact 
solution, we can discuss the collision of black holes 
in the power-law contracting universe just as the 
brane collision [l8| . 

The plan of the rest paper is as follows. We shall 
begin by reviewing the dynamical "black hole solu- 
tion" derived in [2(|. Section Hill involves a detailed 
examination about the properties of matter fields. 
Our main result is contained in Section IIV1 where 
we will elucidate spacetime structures based on lo- 
cal and global perspectives. We draw the confor- 
mal diagram that allows us to pictorially identify 
a black hole embedded in an expanding universe. 
Section [VII summarizes our conclusive results with 
several future outlook. In order to keep the main- 
stream of the text, we relegate some issues to Ap- 
pendixes. 

We shall work in the units of c = K = 1 and 
retain the gravitational constant n 2 = 8ttG. Wc 
pertain to the notation R to denote the circumfer- 
ence radius, so that we use the script throughout 
the paper for the spacetime curvature W 1 vpa V u := 



2V [p V c 



and TZ 



TV, 



II. TIME-DEPENDENT SOLUTIONS 
FROM THE INTERSECTING BRANE 
SYSTEM 

The authors of [2(| have classified the possi- 
ble time-dependent intersecting brane systems in 
M-theory by assuming the metric form, and pre- 
sented some interesting solutions in lower dimen- 
sions by compactification. Among other things, 
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one of the most interesting solutions is a "black 
hole" in the expanding universe. In the case where 
the all branes are at rest, i.e., spacetime is static, 
the 4D reduced solution indeed describes a black 
hole, which is obtained from M2-M2-M5-M5 brane 
system (four brane charges) or from M2-M5-W-KK 
brane system (two brane charges plus a Brinkmann 
wave and a Kaluza-Klein monopole) . The 5D black 
hole solution is similarly derived from M2-M2-M2 
brane system (three brane charges) or from M2- 
M5-W brane system (two branes plus a Brinkmann 
wave). Time-dependent extensions of those lower- 
dimensional solutions are produced from the dy- 
namical intersecting brane systems, in which only 
single brane is time-dependent [2(3]. These 11D so- 
lutions and the procedure of the dimensional re- 
duction are shortly summarized in Appendix |A"1 

After the toroidal compactification of 11D in- 
tersecting brane solutions, we have a 4D "black 
hole solution," which is spherically symmetric, time 
and radial coordinate-dependent (we relegate the 
extension into the case of non-spherical, multiple 
black holes to Appendix |EJ . The 4D metric in the 
Einstein frame is given by 



and 



ds 2 = -Sdt 2 



l(dr 



2 1 ~ 2 dtt 2 ) 



(2.1) 



with 



Qt 
r 

r 



1 



r 

Qs 



-1/2 



> (2-2) 



where dO| = dO 2 



sin Od<j) 2 denotes the line ele- 
ment of a unit round two-sphere. The constants 
Qt and Qs,Qs' ,Qs n are charges of one time- 
dependent brane and three static branes, respec- 
tively. Here and hereafter, the script "T" and 
"S"' are understood to trace their origin to time- 
dependent and static branes. The above metric 
manifests that the conditions of stationarity and 
asymptotic flatness were both relaxed. 

Assuming t/to > and changing to the new time 
slice i defined by 



1.(1) 

to \toJ 



3/4 



with t = -to , (2.3) 



we are able to put the solution (|2.1j) into a more 
suggestive form, 



ds 2 



-SdP + ^t (dr 



2 dCl 2 2 ) , (2.4) 



where 



i + ir 1 1 1 



Qv_ 

r 



Qt 



Qs" 



-1/2 



(2.5) 



to 



1/3 



(2.6) 



When we take the limit of r — > oo, we can find 
that the metric (|2.4[) asymptotically tends to a flat 
FLRW spacetime, 



ds 2 



-dt 2 



(d: 



• 2 + r 2 dn 2 ) 



(2.7) 



Here, the scale factor expands as a cx t 1 / 3 , which 
is the same as the expansion law of the universe 
filled by a stiff matter. Hence we expect that this 
spacetime is asymptotically FLRW universe with 
the equation of state P = p. 

On the other hand, taking the limit r — > with 
t being finite, the time-dependence turns off and 
the metric (|2.4[) reduces to the direct product of 
2D anti-de Sitter (AdS2) space with round sphere, 

ds 2 = -^dt 2 + ^(dr 2 +r 2 dn 2 ) , (2.8) 

where Q := (QtQsQs'Qs") 1 ^ 4 plays the role of 
curvature radii of AdS2 and S 2 . This is a typical 
"throat" geometry of an extreme black hole. To 
take the case of an extreme RN spacetime with the 
Arnowitt-Dcser-Misner (ADM) mass M = Q, the 
metric in the isotropic coordinates reads 



ds 2 



Q 



at 2 



+ 1 + - {dr 2 +r 2 dn 2 ) , (2.9) 



which indeed asymptotes to the spacetime (|2.8p in 
the limit r — > with t kept finite. 1 Thus, we may 
speculate that there might exist a degenerate event 
horizon at r = 0. 

For this reason stated above, it might be tempted 
to regard the present spacetime (|2.1[) as a degener- 
ate black hole with radius Q residing at the center 
r = of the expanding FLRW universe. However, 
there is a prioi no guarantee that the global pic- 
tures of the solution is obtainable simply from the 



It should be emphasized that the AdS2 x S geom- 
etry H2.8H indeed approximates the whole portion of 
the near-horizon geometry-the three dimensional null 
surface-of an extremal RN solution (and more generally, 
of any kinds of an extremal black hole |4(J ) , not restricted 
to the "throat." To the contrary, the spacetime {2.8$ fails 
to describe the near-horizon geometry of the event horizon 
of our time-dependent solution. As proof of the incident, 
we will show in later section that the temperature of a 
dynamical black hole docs not vanish in the present case, 
unlike the extremal RN black hole having a zero Hawking 
temperature. 
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asymptotic considerations, just as we inferred the 
causal structure of the Schwarzschild-AdS space- 
time from those of the Schwarzschild and the AdS 
metrics. Furthermore, it is far from clear to what 
extent the spacetime shares the physical proper- 
ties as its limiting counterparts. In fact, a more 
detailed argument casts doubt on the above naive 
expectation. 

To see this, it is instructive to consider the 
McVittie spacetime [35| . 



d,s 2 = - 



1 - M/2ar 
1 + M/2ar 



dt 2 



2ar v 



r 2 dO|) 



(2.10) 



where a = a(t). The metric (gjjj) with (f2"3j) 
and (|2.6[) looks quite similar to the McVittie space- 
time in appearance. In the limit of r — > oo, the 
McVittic spacetime asymptotes to the flat FLRW 
universe with the scale factor a. If a is set to be 
constant (a = 1), we recover the Schwarzschild 
spacetime with the ADM mass M written in the 
isotropic coordinates. So one might deduce that 
the McVittie spacetime (|2.10[) describes a black 
hole immersed in the FLRW universe. 

However, it turns out that this exemplifies that 
our naive estimate ceases to be true. 2 When the 
scale factor obeys the power-law form a oc t p (p > 
0), a curvature singularity appears at r = M/(2a). 
As Nolan pointed out [3g, the spacetime events 
described by r = Mj (2a) in part consist of a shell- 
crossing spacelikc singularity lying a future of a big 
bang singularity a = 0. The surface r — Mj (2a) 
fails to describe a (regular) horizon of a black hole 
in FLRW universe. Inferring from the "monopolc 
term" m/r, docs this spacetime instead describe a 
point mass (singularity) at r = embedded in the 
expanding universe? The answer is NO. It turns 
out that r = corresponds to infinity, rather than 
the locus of a point particle in the universe. Besides 
that, according to the quasi-local definition of hori- 
zon, the McVittic spacetime may serve as a "white 
hole" in an expanding FLRW universe. As a good 
lesson of above, we are required to take special care 
to conclude what the present spacetime describes. 

In this paper, we study the above space- 
time (|2.ip more thoroughly [we are working mainly 
in Eq. (|2.1j) rather than Eq. (|2.4p . because the for- 
mer coordinates cover wider range than the latter] . 
We assume to > 0, viz, the background universe is 



2 The case a = e Ht with constant H is exceptional for which 
the metric describes an SdS spacetime. But we are reluc- 
tant to refer to this as a black hole in expanding universe, 
since the metric is rewritten in a static form thanks to 
Birkhoff 's theorem. 



expanding. For simplicity and definitcness of our 
argument, we will specialize to the case in which all 
charges are equal, i.e., Qt = Qs = Qs' = Qs" = 
Q (> 0). 3 To be specific, we will analyze the space- 
time metric 

ds 2 = -Edt 2 + (dr 2 +r 2 dnl) , (2.11) 

whose component S, Eq. (|2.2[) . is simplified to 



5 - (H T H^) 



with 



t±T — 1 , 



He 



1 + - 

r 



(2.12) 



(2.13) 



A more general background with distinct charges 
are yet to be investigated. The result for the 5D 
solution will be given in Appendix [Cj 



III. MATTER FIELDS AND THEIR 
PROPERTIES 

It is a good starting point to draw our attention 
toward the matter fields. Since we know explicitly 
the 4D metric components, we can read off the total 
energy-momentum tensor of matter fluid(s) from 
the 4D Einstein equations, 



k 2 T u 



(3-1) 



where = 1Z^ V — {lZ/2)g^ u is the Einstein ten- 
sor. What kind of matter fluids we expect? There 
may appear at lest two fluid components: one is 
a scalar field and the other is a U(l) gauge field. 
This is because we compactify seven spaces and 
we have originally 4- form field in 11D supergravity 
theory. The torus compactification gives a set of 
scalar fields and the 4-form field behaves as a U(l) 
gauge field in 4D. In our solution, wc assume four 
branes, which give rise to four U(l) gauge fields. 

As shown in Appendix [XJ we can derive the fol- 
lowing effective 4D action from 11D supergravity 
solution via compactification, 



S 



2 k 



-ft--(V<I>) 2 



16tt ^ 



(3.2) 



' ! If Qt is different from other three same charges Qs, the 
present result still holds. It is because such a difference 
amounts to the trivial conformal change 

ds 2 = (Qt/Qs) 1/2 [-H*diJ +S7 1 (dr 2 +dfi|)] 

with simple parameter redefinitions: H„ = [(i*/t*rj + 
?s/r)(l + Qs/r) 3 ]- 1 ' 2 , U = (Q T /Qs)- 1/2 t, and i.„ = 



(Q T /Qs) 1/2 to- 
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where $, and X A {A = T,S,S',S") are a 

scalar field, four U(l) fields, and coupling con- 
stants, respectively. 

The above action yields the following set of basic 
equations, 



7 » ^ A ^ F (A)^ =0t 



where 



T-'(cm) 



47T ^ 



F (A)p(A)p 
fJ-P v 



(3.3) 
0, (3.4) 

(3.5) 
(3.6) 
(3.7) 



For the present case with all the same charges, two 
different coupling constants appear. 

A simple calculation shows that our spacetimc 
metric (|2.11j) satisfies the above basic equations 
()3.3j) . (|3.4j) . and (|3.5|) . provided the dilaton profile 



(3.8) 



and four electric gauge-fields 

(T) _ \ / 2^Q 

^oi - - 



r 2 iJ 2 ' 



kF, 



(S) 



kF (s') _ F (s") 
K_r 01 — K_r 01 



V2nQ 



(3.9) 



with the coupling constants 
At = v6, As = Ag' = As« 



V6/3. (3.10) 



The U(l) fields are expressed in terms of the elec- 



trostatic potentials = V ^Av — V U A^ as, 



kA, 



IT) 



(3.11) 



as 



tis 



where we have tuned A^ to assure Aq S ^ 
r — > oo using a gauge freedom. Therefore the 
present spacetime (|2.11j) is the exact solution of 
the Einstein-Maxwell-dilaton system Q3.2p . 

One may verify that Q is the (electric) charge 
satisfying 

7& = ih j/ A ^ F ^ )AS "^ ( 3 - 12 ) 

where S is a round sphere surrounding the source. 
This expression is obtainable by the first integral 
of Eq. (|3~5)) . 

Note that one can also find magnetically charged 
solution instead of (|3.9[) . However, this can be re- 
alized by a duality transformation 



(3.13) 



which is a symmetry involved in the action 
Henceforth, we will make our attention only to the 
electrically charged case. This restriction does not 
affect the global spacetime picture. 



A. Energy density and pressure 

Using our solution (|2.1[) . we can evaluate the 
components of the energy- momentum tensors, i.e., 
the energy density and pressures for each field ( the 

dilaton $ and U(l) fields F^ } ). They are given by 



~ <J>' 2 



p(<&) = p ($) 



1 



5"1 $2 _ s$ /2 



2 

(cm) _ _p(cm) _ p(cm) _ p(cm) 



16k 


2 


3 




16k 


2 


1 


> 


8^ 




Q 2 



tl \R 



T 



3/2 g 2 (g T -g s ) 2 " 

n 3/2 Q 2 (ff T ~ff s ) 2 " 



> 0. 



1 



5 Xt "*(F ( P) 2 + 3e AsK *(F ( f } ) 2 



4 K 2 r 4 



H T 



3/2 



-1/2-1 



> 0. 



(3.14) 
(3.15) 

(3.16) 



The time t — to is special at which the energy 3/(16^), and soon after it becomes gradually in- 
density of the scalar field is uniform K 2 p^(to) = 



G 



homogeneous. We shall refer to to as the fiducial 
time. As the universe expands, the energy den- 
sity of the scalar field at infinity behaves p'®' oc 
t~ z l 2 oc i~ 2 oc a~ 6 , as expected for the FLRW uni- 
verse with stiff matter or a masslcss scalar field. 
The energy density of the U(l) fields evaluated at 
the fiducial time is «V em) (*o) = Q 2 /(rH s ) 4 = 
Q 2 /(r + Q) 4 , which is the same as that of the 
extreme RN spacetime, and decreases in time as 
K 2 p( em )(t, r) oc t~ x l 2 oc t~ 2 / 3 oc a" 2 near infinity 
as in the same manner for the FLRW universe. 

We can also find the energy flux J M := T 0ai whose 
spatial component J- is given by 



k 2 T 



k 2 T\ 



K 2 T (*)6 



3Q (H T - H s ) 
8t r 2 H 2 H s 



(3.17) 



where a hat denotes the tetrad component. Only 
the scalar field contributes to the energy flux, since 
there exist no magnetic fields, i.e., no Poynting 
flux. One can find that the T becomes negative 
for t > to j implying that the scalar field energy is 
falling toward the black hole. However, as it turns 
out later, the flux never gets into the black hole. 
This may be attributed to that the repulsive force 
caused by the U(l) fields becomes strong near the 
horizon and finely balances the attractive gravita- 
tional force of the dilaton field. 

It is worthwhile here to discuss the issue of the 
energy conditions. Many candidates of "black hole 
solutions" in expanding universe found in the liter- 
ature (McVittie's solution [35|], Sultana-Dyer solu- 
tion [34| . etc) do not respect energy conditions in 
the whole of spacetime. Now the present system- 
U(l) gauge fields coupled to dilaton-apparently 
satisfies the energy condition and hence it provides 
us a nontrivial example. We may verify this explic- 
itly as follows. Inspecting Eqs. (|3.14[) - p.l6[) . we 



notice p^ = P r 



> P^ ] and p(° m ) 



-P ( 



> 



P, 



(cm) 



from which we obtain 



p =p m + p {cm) = P r w - P r (om) 
>P r w + P 7 < cm) =P r , 



p =p(*> + p^ > Pf ) + P 9 (cm> = Pe , (3.18) 
and 

p+P r = {pW + P r W ) + ( P {Cm) + P r (Cm) ) > , 

p + P e = (p<*> + Pf >) + (p (cm) + P e (cm) ) > . 

(3.19) 

These equations mean that p > |Pj| (i = 1, 2, 3) is 
satisfied anywhere. 

The energy flux = T 0/i satisfies 

2 , x-2 



J^J M = —pr + T 



-(p( om )) 2 -2p(° m )pW-(pf ) ) 2 <0 



where we have used the relation (p^) 2 — T 2 = 
(Pg*-*) 2 at the second equality. Hence, the energy 
current J M is a timelike vector everywhere. It then 
follows that the spacetime (|2.1|) satisfies the domi- 
nant energy condition [p > \Pi\ (t = 1, 2, 3) and J M 
is non-spacclike] . 



B. Misner-Sharp energy 

Another useful quantity to characterize matter 
fields is the quasilocal energy, which is defined on 
the closed two-surfaces. If the spacetime has an 
SO(3)-symmetry, we are able to give a physically 
satisfactory quasilocal energy introduced by Misncr 
and Sharp [|l|. 

The utility of spherical symmetry lies in the fact 
that we can covariantly employ the circumference 
radius, 



R(t,r) 



IrlS- 1 / 2 



Irl(PTPi) 1 / 4 , (3.21) 



(3.20) 



in terms of which the area of metric sphere is given 
by 4ttR 2 . R is a geometrical quantity and has an 
invariant meaning. It is occasionally of great ad- 
vantage, instead of the comoving coordinate r, to 
make use of R. Using the circumference radius, the 
Misncr- Sharp quasilocal energy is defined by [4l|, 

A-kR 

m(t, r) := — [1 - g" v (W„R)(W v R)) . (3.22) 

This quantity is a useful local measure to demon- 
strate geometric properties of spacetime j42Tj4^ |. 
The Misner-Sharp energy represents a mass energy 
contained inside the surface of radius R. Once 
the compact surface is specified, the Misner-Sharp 
mass is given without any ambiguity. Such a quasi- 
localization is possible because of spherical symme- 
try, in which no gravitational wave exists. By def- 
inition, it is characterized by geometric structure 
(metric components and its first derivatives) and 
does not require the premise of asymptotic struc- 
ture of spacetime. So, it is considerably advanta- 
geous for the analysis of local spacetime structure. 

Physical interpretation of the Misner-Sharp en- 
ergy has been further backed by various desirable 
properties: it satisfies the first law of thermody- 
namics [ill , and it shows the properties of positiv- 
ity and monotonicity under the dominant energy 
conditions, and it reduces to ADM mass in the 
asymptotically flat spatial infinity. Reference [43j 
reinterpreted it by the integral of a locally con- 
served energy current coming from the symplectic 
structure of spherical symmetry. One may find the 
superiority of the use of Misner-Sharp energy in the 
next section. 

The present metric (|2.11l) with (|2.13[) gives rise 

to 
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_ Tt\r\[r 2 H T Hl + l6t 2 H 2 H 2 - (H s t + 3H T t a ) 



..,., L . , , 

4^(^^5)1/4 v ; 

Physical meaning of each term in this equation is best understood at the fiducial time t = to, at which 
the Misner-Sharp mass is expressed as 



*»mfo, R) = tMjW + + 4vrQ = «' ^ 



R 3 P m {t )+41T J i? 2 di?p(° m )(i ,i?)^ +47TQ. 



(3.24) 



The first term corresponds to the energy of the 
scalar field, and the second and the last terms to 
the U(l) energies outside of and inside of the black 
hole, respectively. 

If we set Q = 0, the second and third terms 
in Eq. (|3.23[) are combined to cancel, and m is 
expressed by the coordinates (12. 4[) as 

ar 3 ( da\ 2 4tt 3 (4>) 
m= 2GUJ = Y (-)P ( >, (3-25) 

as expected for the background FLRW universe. 
This also justifies the first term in (|3.23|) to be the 
contribution of the scalar field. 



IV. SPACETIME STRUCTURE 

Let us now turn to the main task of revealing 
spacetime structure of the solution (|2. 1 1|) . In order 
to address this issue, it is of significance to discuss 
the followings: 

• Singularity 

• Trapped surface 

• Event horizon 

• Asymptotic structure 

The first two topics are associated with the local 
character of spacetime, whereas the last two re- 
quire global considerations. These are elementary 
issues to be explored in order to characterize the 
spacetime. 

We wish to show that our metric (|2.11[) describes 
a black hole in the FLRW universe. To this end, 
we need to establish in the first place that the 
far outside region from the central inhomogeneous 
domain behaves non-pathologically. If the space- 
time admits a naked singularity in the asymptotic 
region-other than the initial big bang singularity- 
that are not covered by the event horizon, the so- 
lution would not gain a popularity as a black hole. 



According to a series of theorems due to Pen- 
rose and Hawking [22| . the appearance of space- 
time singularity is closely associated to the pres- 
ence of trapped surfaces. Thus, the examination of 
trapping property may provide us useful informa- 
tion of local spacetime geometry. In particular, we 
have two competing effects due to the black hole 
and the expanding cosmology: the former tends to 
focus light rays back into the hole while the latter 
tends to spread it out to infinity. 

At first sight, one might expect to gain only a 
limited perception from the local point of view, 
even though the curvature singularity and local 
horizons are indeed of importance. Nevertheless, 
under physically reasonable circumstances, the ex- 
istence of the local horizon implies that the event 
horizon lies outside it (22j. To be more precise, all 
trapped surfaces are contained within black holes 
under the null energy condition provided the space- 
time is asymptotically flat with some additional 
technical assumptions (see Proposition 9.2.8 in [22j 
for the proof). So, this criterion is of use for our 
study. 

In order to define a black hole as a "region of 
no escape," the spacetime must allow null infinity 
as an idealization of an observer sufficiently far in 
the distance from some central region. Thus the 
asymptotic analysis should also be put on an em- 
phasis. In particular, we will take a close look at 
null geodesic motions, since the null rays play a 
privileged role in the black hole geometry. Above 
listed considerations are sufficient to provide us 
with insight into the global pictures of our dynam- 
ical spacetime (|2.4[) . 

In the remainder of this paper, we shall simplify 
our notations by using dimensionless variables t := 
t/to and r := r/Q. The metric is also rewritten 
into dimensionslcss form ds 2 = Q~ 2 ds 2 as 

d§l = -T 2 § dP + S" 1 (dr 2 + f 2 dn 2 2 ) , (4.1) 

with 

r .= *° S. = f l (42) 

Q ' [(l + if)(l + f) 3 ]V2 ■ { ■ ' 
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In this form, the metric involves only one dimen- 
sionless parameter r = to/Q, the physical mean- 
ing of which will be revealed below. We will affix 
"tilde" to denote dimensionless quantities in what 
follows. 



7? 7?^ 



These curvature invariants diverge when Ht = 
and Hs = 0, that is, at 

t = t s (f) := -l/f, and f = -1 . (4.5) 

At these spacetimc points, the circumference radius 
R, Eq. p.21j) . vanishes, i.e., they are central shell- 
focusing singularities. Thus, around infinity is free 
from singularities and is well-behaved. 

It deserves to observe that t = surface, where 
the scale factor a(t) appearing in the metric (|2.4[) 
vanishes, is not singular at all since the curvature 
invariants remain finite therein. It follows that the 
big bang singularity t = is smoothed out due to a 
nonvanishing Maxwell charge Q (> 0). Hence, one 
has also to consider the t < region in the coor- 
dinates (|2.11[) . In addition, we find that the f = 
surface is neither singular, thereby we may extend 
the spacetime across the f = surface to r < 0. 
Since the allowed region is where HtH^ > is 
satisfied, we shall focus attention to the coordinate 
domain 

t>t 3 (f), r>-l, (4.6) 

in the subsequent analysis. Another permitted re- 
gion t > t s and f < — 1 is not our immediate in- 
terest here, since it turns out to be causally dis- 
connected to the outside region, as we shall show 
below. Possible allowed coordinate ranges are de- 
picted in Figure [T] 

Since our spacetime is spherically symmetric, 
electromagnetic and gravitational fields do not ra- 
diate. Thereby, it is more advantageous to con- 
centrate on their "Coulomb components." For this 
purpose, let us introduce the Newman-Penrosc null 



A. Singularity 

The scalar curvature and the Kretschmann in- 
variant scalar are given by 



r 




FIG. 1: Allowed coordinate ranges. The grey zone 
denotes the forbidden region, and the dashed curves 
correspond to curvature singularities. 



tetrads by 

Z M dx" = ^(-rdi + H-Mf), 

n^dx^ = y|(-Td£- S _1 df), (4.7) 
v 

m^dx^ = ^_ (dO + i sin 6&4>) . 

with fhft being a complex conjugate of m^. They 
satisfy the orthogonality conditions l^n^ = — 1 = 
—m^rh^ and l^l^ = n^n^ = m^m^ = fhP'rh i i = 0. 
Since t is a timelike coordinate everywhere, Z M and 
arc both future-directed null vector orthogonal 
to metric spheres. 

The only nonvanishing Maxwell and Weyl scalar 

are their "Coulomb part," (/)[ A) := -^F^ (l^n" + 
rh^m 11 ) and ^ := —C^ vpa l^m v fhPri'', both of 
which are invariant under the tetrad transforma- 



p. 



3 

w 

T 2 rrlO 



1_ /Ht 
t 2 \H S 



3/2 



1 (H T ~ Hs) 1 



f 4 {H T H 7 s y/ 2 

T2 



(4.3) 



15H T H^ + 6 ^ H T H* S (7H T + 1QH T H S - ffj) 
{9(317J| + 2H S + 159)H T + 12(7H% + 9H S + 96)H T H S 
6(15iff - 1267Js + 143)7J|7J| - 12(H S - 1)(H S + 15)7J T i/| + 7l(H s - l) 2 H s } 



(4.4) 
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tions due to the type D character. It is readily 
found that 



j(T) 



V2nQ? 2 H 2 ; 



b (s) 



V2nQ? 2 H 2 s 



(4.8) 



and 

*2 = 



i — i r Ts — i . 



i 2 + 1trH% 



6Q 2 



(4.9) 



The loci of singularities at which these quantities 
diverge are the same as the positions of the above 
singularities. One may also recognize that at the 
fiducial time t = 1, above curvature invariants are 
the same as the extremal RN solution, as expected. 

Let us next look into the causal structure of sin- 
gularities. Since the Misner-Sharp mass p.23[) be- 
comes negative as approaching these singularities 
[the third term in Eq. (|3.23j) begins to give a dom- 
inant contribution], we speculate from our rule of 
thumb that these singularities are both contained 
in the untrapped region and possess the timclikc 
structure. 

To see this more rigorously, let us consider radial 
null geodesies 



df 
dl 



±rS. 



(4.10) 



in the neighborhood of these singularities. Here, 
the upper (lower) sign corresponds to the geodesies 
along the direction l>* {n u ), which we shall refer to 
as outgoing (ingoing). If we can find an infinite 
number of null geodesies that emanate from and 
terminate into the singularity, the singularity turns 
out to be timelike. If there exists a unique geodesic 
terminating into (emanating from) the singularity 
and an infinite number of geodesies emanating from 
(terminating into) it, the singularity has an ingo- 
ing (outgoing) null structure. Whereas, if there 
exist an infinite number of ingoing and outgoing 
geodesies terminating into or emanating from the 
singularity, we can conclude that it is spacclike. 



J 



We begin by the analysis of the singularity at r = 
— 1. In the vicinity of the singularity f — — 1, we 
suppose that the null geodesies have the following 
asymptotic solution, 



1 = Ci \t-tx 



(4.11) 



where Ci, t\ and p are constants. Ci is taken to 
be positive since we are concerned with the region 
f > — 1. ti can be regarded as the arrival time 
of the ingoing null geodesies at the singularity, or 
the departure time of the outgoing null geodesies 
from the singularity. Consider first the t > t± case, 
that is, the geodesies emanating from the singular- 
ity. Substituting the assumed form (|4.11|) into Eq. 
(14.1 Op , we find that only the outgoing null geodesies 
have the solution, for which 



5/2 



-1/2 



(4.12) 



This reveals that for t > t\ there exist radial null 
geodesies that departed the singularity at t = t\. 
Following the identical procedure, it can be shown 
that for t < ti only the ingoing null geodesies have 
the solution for which C\ and p are given by (|4,12[) . 
This means that there exist radial null geodesies 
that will get to the singularity at t = t±. There- 
fore, we establish that the singularity f = — 1 is 
truly timelike, that is, a locally naked singularity, 
since we have only one set of outgoing and ingo- 
ing null solutions parametrized by their arrival or 
emanating time. It is also obvious from Eq. (|4.12l) 
that the singularity f = — 1 exists only for t < 1. 

In an analogous fashion, assume the asymptotic 
form of the null geodesies near the singularity r = 
r s (= -1/t) as 



l=C 2 \t-t 2 \ 



(4.13) 



where C 2 , t 2 and q are constants. In the present 
case, C 2 takes positive (negative) value for t 2 > 
(h < 1). 



1 



Plugging this into Eq. (|4.10[) . we obtain q = 2/3 and 



cr = 


3 

— r 

2 


[*a(*2 ~ 


-i) 3 ] 


-1/2 

) 


for outgoing null with t > i 2 , 


(4.14a) 


c 2 f' 2 = 


3 

— T 

2 




hf\ 


-1/2 

) 


for ingoing null with t > t 2 , 


(4.14b) 


cr - 


3 

— T 

2 




-if] 


-1/2 

) 


for ingoing null with t < t 2 , 


(4.14c) 


c 2 f' 2 = 


3 

2 T 


m- 




-1/2 

) 


for outgoing null with t <t 2 . 


(4.14d) 
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Eqs. (|4.14b|) and ()4.14dj) indicate that there exist 
null geodesies that originate from and come to the 
singularity at any time £2 < 0. Thus, the singular- 
ity t = —1/f (< 0) appearing in the f > region 
is also timelike. Whereas, Eqs. (|4.14aj) and (|4.14cj) 
show that the timelike singularity r = — 1/t (< 0) 
occurs only for t > 1. 

Although we have not dwelt on the existence 
proof of the asymptotic solutions [Eqs (|4.11[) with 
(|4TT2|) and (|4~T3)) with (|4~T4]) ] of the radial null 
geodesies, the proof can be obtained via the con- 
traction mapping method following the argument 
in e.g., the appendix of [46|. It is also noted 
that it is sufficient to focus attention on the ra- 
dial null geodesies in the vicinity of central sin- 
gularities: causal geodesies excluding radial null 
geodesies will fail to emanate from singularity if 
radial null geodesies do not arise (see also the ap- 
pendix in |46[). 



B. Trapping horizons 

Our primary concern in this article is to reveal 
that the spacetime describes a black hole. How- 
ever, this is not only technically but also concep- 
tually difficult. Since a black hole is by definition 
a "region of no escape," the locus of event horizon 
as its boundary has a teleological meaning. Unless 
we know the entire future of our universe, we are 
unable to determine whether a "black hole candi- 
date" is qualified black hole. The concept of 
the black hole is considerably messy in practical 
point of view. 

In physically acceptable situations, however, we 
can rely on the the notion of trapped surface, first 
introduced by Penrose [!?} . The concept of trapped 
surface is inherently local, such a difficulty does not 
arise. Imagine a massive star undergoing a gravita- 
tional collapse to form a black hole. There appears 
a region for which even "outgoing" null rays are 
dragged back due to strong gravity and have a neg- 
ative expansion. For each time slice, this defines 
an apparent horizon p2j as an outermost bound- 
ary of the trapped region in the asymptotically flat 
spacctimcs. Hayward generalized these quasilocal 
concepts to define a class of trapping horizons [48j . 
One strength of the use of trapping horizons is just 
to encompass various types of horizons associated 
not only with black holes but also with white holes 
and cosmological ones. As commented at the be- 
ginning of this section, the underlying aim of this 
direction is to gain useful guide for event horizon 
from these local analysis. 

Let us consider a compact spacelike orientable 
surface S. We take S as a metric sphere, respect- 
ing an SO(3)-symmetry of background universe. 
It then follows that the Newman-Penrose tetrads 
(P, n**) defined in Eq. (|4. 7[) are normal to S (i.e., 



they are radial) and future-directed null vectors. 
Due to the spherical symmetry, they are shear-free 
and rotation-free. Define the associated null ex- 
pansions 9± by 

9+ := 2m ( ' 1 m I/) V^ i „ 9- := 2m (At m I/) V A1 n„. 

(4.15) 

In the coordinates (|4.ip . they are expressed as 
~ _ rggCgrgf)Va ± t(3H t + H s t) 

± 2V2Tf{H^H 7 s y/ i ■ y • J 

Note that the signs of 9± have an invariant mean- 
ing, but each value of 9± is not a universal quan- 
tity due to the "class III" tetrad rotations. An 
invariant combination is their product 9 + 9- = 
-2i?- 2 (V Al i?)(V' I i?). 

Expansions 9± characterize the extent to which 
the light rays are diverging or converging, or equiv- 
alently the rate at which the area of metric sphere 
is increasing or decreasing in the null directions. In 
terms of null expansions 9± , a metric sphere is said 
to be trapped (untrapped) if # + 6L > (9+9- < 0), 
and marginal when 9+9- — 0. A marginal surface 
is said to be future (past) if 9 + = (6- = 0). A fu- 
ture marginal surface is further classified into outer 
(inner) if n^V M #+ < (n''V ^9+ > 0) and similarly 
a past marginal surface is called outer (inner) if 
l^V ^9- > (l^y^O- < 0). A trapping horizon 
is the closure of a hypersurface foliated by future 
or past, and outer or inner marginal surfaces [48j . 
In terms of the Misncr-Sharp energy, R < 2Gm 
(R > 2Gm) defines the trapped (untrapped) region 
and the marginal surface (9+9- = 0) is positioned 
at the "Schwarzschild radius" R = 2Gm. A "nor- 
mal" spacetime region as occurred in the flat space 
is composed of untrapped surfaces on which outgo- 
ing rays have positive expansions while the ingoing 
rays have negative expansions. 

Among these classes of trapping horizons, the 
future-outer trapping horizons turn out to be most 
relevant in the context of black holes. The future 
outer trapping horizon properly captures the intu- 
itive idea that the ingoing null rays are converging 
with the outgoing null ray being instantaneously 
parallel on the horizon, diverging outside and con- 
verging inside. Inner trapping horizons are associ- 
ated with cosmological horizons, and interior hori- 
zons of black holes. The past trapping horizons 
arise when discussing white holes and cosmologi- 
cal ones. Since the concept of trapping horizons 
is sufficiently general, it is considerably useful for 
the analysis of black holes especially in the non- 
asymptotically flat spacetimes. 

The expansions arc intimately associated to the 
variation of the Misner- Sharp energy through the 
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first law [42|, 

fV^m = 2nR 3 {T^n u 9^ 
n^V M m = 2tt R 3 (T l ,J> J -n v 6- 



- T^n v 9+) . 



(4.17) 



The present system satisfies the dominant energy 
condition, which implies T^l u > 0, T ilv n ii n v > 
and T^J^n" > 0. Hence, Eq. (|4~T7|) establishes 
that m is not decreasing (not increasing) along the 
^-direction (n^-direction) in the untrapped region 
of 9+ > and 6- < 0. This illustrates that the 
Misncr-Sharp energy is a monotonically increasing 



function toward outwards in the ordinary region 
where 9+ > and < are satisfied, which ac- 
cords with our intuition. 

Let us now investigate the properties of trapping 
horizons more closely. From Eq. (|4.16|) . trapping 
horizons occur at 

fH s {H T H 3 s ) 1/2 ± t(3H t + H s t) = 0, (4.18) 



where the upper (lower) sign corresponds to 9 + = 
(0_ = 0). Noticing that rH s > 0, Eq. (|4TTg|) is 
solved as t — fc^j (r), where 



J 



f 2 



2r 2 ( J H s + 3) i 



H 5 s -6t 2 (H s 



3)r- 3 T #1 



^1 



4t 2 (# s + 3)f 



r-3 



(4.19) 



r 



Here, 9± = holds at t = t^(^). In the pursuing 
subsection, we shall separately analyze the r > 
and f < cases corresponding to outside and inside 
the black-hole event horizon. 

Before going into the detailed argument, wc 
pause for a moment to discuss the behavior of 
R(t,f). We have taken P (and correspondingly 
n M ) in such a way that r increases (decreases) along 
P (n* 4 ). From 



This means that R increases as f grows for t > t c , 
deserving to be called "outgoing." However, this 
is no longer true for t < t c (f). This means that 
there exists a maximum value of R(t, r) for given 
time (see Figure [2]) . It should be also remarked 
that the invariant scalar curvatures are all finite at 
this surface t = t c (r), this is not the shell-crossing 
singularity. 



R ,t 



H 3 
H 3 



1/4 



R f 



\f\[3H T + H s t] 



4r(H*H s ) 



1/4 ' 

(4.20) 



one finds that R.f > (R t f < 0) holds for t > t c (f) 
(t < t c (r)), where 



1. Trapping horizons in the region off > 



4f+ 1 



(4.21) 



Let us begin by the case of r > 0. After simple 
calculations, one obtains 



r(+) 



TH 



t c (r) — i TH 
?(-) 



3)2 \ h3 s\/ H I + ^r- 3 (H s +3) + Hl)>0, 



t 



TH 



ts(r) 



2t(H s + 3) 2 
f 2 

2t(H s + 3) 2 



H s\^ H s + 4r 2 f-3(7? s + 3) - H% > 



1 



H s [Hi- JH* + At^- 3 {H s + 3) 



(4.22) 
(4.23) 

>0, 
(4.24) 



for f > 0. Figure [3] shows the typical curves of denotes a past trapped region in which even in- 
trapping horizons . The region < * < *th going light rays are diver S in g due to the cosmic 
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R 



3.0 






15 


t=2y 










2.0 






15 













1/2 



FIG. 2: Typical behaviors of R = [(tf+l)(f +1) 3 ] 1/4 as 
a function of f with fixed time. For fixed t < 0, the cir- 
cumference radius R ceases to increase monotonically 
with f , but has a maximum at t = t c (r). For t > 0, R 
turns to be monotonic in f. 




FIG. 3: Typical curves of trapping horizons £th(?) f° r 
T < r cr it outside the event horizon f > 0. The plot 
is t = 1. tj^j first occurs at some time (in this case 
t = f = 1) and bifurcates into two branches. The black 
dashed line and the grey dotted line denote t 3 (f) and 
t c (r), respectively. 



trapping horizon occurred in the r > 0, t > 
region is always of outer-type. This may be as- 
cribed to the cosmological origin since a past-outer 
trapping horizon develops when the background 
FLRW universe is filled with a stiff matter. More 
stringent bound for the condition F > is numer- 
ically found to be t < r cr i t ~ 5.444. In the case of 
t > T cr it i a part of the past trapping horizon near 
f — becomes inner rather than outer. 
Similarly, we obtain 



F 



4f 4 (7J5F|)i/2 



t-t ( -' 

TH 



(4.28) 



along the trapping horizon with 6 + = 0. From Eq. 
gJH , we find 

F| t - =t -(-) (?) > flf - (4r~ 2 + 8r + 1)H^ 

= 12H 2 S (3 + H s y 2 > 0. (4.29) 

Hence, the trapping horizon is necessarily of 
future-inner type. The appearance of untrapped 
region 9 + < and > is due to the repulsive 
nature of the timelike naked singularity. 

From the general argument, the outer (in- 
ner) trapping horizons must be non-timelike (non- 
spacelike) under the null energy condition [13, EH . 
So the above analysis asserts that the trapping 
horizon never becomes timelike for r < T cr ; t , 
while cannot be spacelike anywhere. 

Here, it is worthwhile to remark the behavior of 
trapping horizons in various asymptotic limits. In 

the r -> limit, Eq. (|4.19[) yields that ft^H is nnite 
and given by 



expansion. 



Next, let us delve deeper into the type of trap- 
ping horizon. It is straightforward to find that 



F 



t -(+) 

' TH 



(4.25) 



rt {+) - 



1 + VI + 4t 2 
2^2 



1 - v/l + 4r 2 
2^2 



> 0, 



<0. 



(4.30) 
(4.31) 



as f — > 0. Using this equation, one can find that the 
circumference radii of the trapping horizons £^ 
respectively approach to some constants R± as r — > 
0, where 



along the trapping horizon with (9_ = 0. Here, we 
have introduced a function 



F := Hi 



(4f 2 + 8f + 



T • 



(4.26) 



the sign of which controls the type of trapping 
horizon. For f > 0, the inequality f + 1 < 



V4f 2 



8f + 1 holds, from which wc obtain 



F < -t(H T + l/r)(4f 2 



1) < 0. 



(4.27) 



Hence, for the trapping horizon with > 0, 
i M V„0_ < is satisfied. That is to say, the past 



R± 



Vl + 4r 2 ± 1 
2t 



1/2 



(4.32) 



These surfaces correspond to the infinite redshift 
(t oo) and blucshift (t — oo) surfaces with 
respect to an asymptotic observer. We will see in 
the next section that these surfaces represent the 
black hole and white hole horizons. Notice that R + 
(and respectively i?_) is a monotonically decreas- 
ing (increasing) function of r and they behave as 
-R+ — > oo and i?_ — > in the limit r — > 0, while 
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they asymptotically tend to unity as t — > oo. Ac- 
cording to Eq. (|4.32[) . Q and r are expressed in 
terms of R± as 



Q — y/ R+R-, 



R+R- 
Rl-R 2 



(4.33) 



Hence, we find that the charge Q sets the geomet- 
rical mean of horizon radii and their relative ratio 
is encoded in the parameter r. 

The physical meaning of r is found by evaluating 
the energy densities of the dilaton field and U(l) 
field at the horizon R + as follows: Those densities 
are given by 



2 ($) I 
« P { '\R+ 



mlR% 



*V em) U + = 



l + 3i?S 



Then, the ratio is found to be 

p( cm ) _ 2t 2 (1 + 3/^_) 



(4.34) 
(4.35) 
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3R% 



(4.36) 



From the expression of we find that r has a 
one-to-one correspondence to p( em )/p(*) \ R+ and is 
given by 









3 p (o " 


i) 


8 p(*> 






p(*> 







/.(■ 


p(om) \ 

«+ 7 


for 


p (cm) 
P W 


»2, 


for 


p<™> 


- 2. 



as (f 2 /4r) 2 , whence R(t^(f),f) -t r 3 /Ar 
— 1/f = t s (f) as f 



DC'. 



?(-) 



fl(«y(f),f) 



(4.37) 

r is a monotonic function of the ratio, and it 
vanishes when the ratio approaches 2. The ratio 
(p( em )/p( 4> ))| fl+ , which must be larger than 2 for 
t 2 > 0, corresponds to r by a one-to-one mapping. 
Hence r is related to the ratio of two densities at 
the horizon. 

In the limit f — > 0, one finds that 0+ and 6L 
both vanish, implying that the f = surface be- 
comes degenerate into an ingoing and outgoing null 
structure. 

On the other hand, taking f — > oo, 2™ diverges 



While wc have tip H ' — ► — l/r = t s (?"J as r — s- oo. 
But this does not mean that the trapping hori- 
zon plunges into the singularity t s as f — > oo. 
One can verify by taking higher order terms into 
account that i^ii tends to have a constant radius 
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FIG. 4: Typical curves of trapping horizons t TH -(f) for 
T < r cr it inside the event horizon f < 0. The plot is 
for t — 1, for which the trapping horizon has a 
maximum at f ~ —0.1080, spacelike for —0.1299 < f < 
and merge with at fo — —0.2679. 



2. Trapping horizons in the region of r < 



For negative values of r, t^H > 



4"h holds > 



contrast to the f > case. Two trapping horizons 
t^j develop only in the region f o < f < and they 
coincide at the point f = fo, where the square-root 
of Eq. (|4.19[) vanishes. Namely, fo satisfies 



/(f ) := H 4 s (f ) +4T 2 (H s (f ) +3)t=, 



v 3 = o. 

(4.38) 



yfr in this limit. 



Since 7(-l) = -12r 2 < and I (-1/4) = 81 > 0, 
— 1 < ?~o < —1/4 is concluded. We can also find 
that the trapping horizon negatively diverges 
at r = —1/4. Figure 2] plots typical curves of trap- 
ping horizons occurred in f < 0. 

Equations (|4.25l) and (|4.28[) continue to be true 
for f < 0. Wc can find numerically that at t^(f) 
F < holds around f ~ implying that it is space- 
like. The future trapping horizon £pjj(r) an d other 
portion of past trapping horizon ti^(r) become 
timelike. In the limit f — > 0, the trapping hori- 
zons i™(r) have constant circumference radii R± 
as outside. 



3. Constant R surfaces 

The trapping horizons occur at (V AI i?)(V M i?) = 
0, i.e., the surfaces of R = constant become null. 
Then the contours of circumference radius will help 
us to recognize the positions of trapping horizons 
in terms of circumference radius. 
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Solving Eq. (|3.2ip with respect to t, we obtain 



- 1 
t = - 

r 



(1+f) 3 



1 



(4.39) 



Taking the derivative of this equation with fixing R, 



we find the relation between dt and dr. Inserting 
this to the metric (|4.1[) , the line element of R = 
constant surface is given by 



J 



where 



dS 2 



t 2 (1 + 4f) 5 
~" 2 R 2 {l + r) 



Rl = R\t^{f),f) 



Rl = R 2 {t { -^{f),f) 



R — Rf) [R + Rf) [R — R%) [R + R%) dr 



(1+f) 4 
2r|l + 4f| 

(1 + ^) 4 
2r(l+4f) 



, , 4r2(l + 4f) 

1 + V 1+ (i + f)* 



/ 4r2(l + 4f) 

" 1+ v 1+ (i+f) 4 



(4.40) 

(4.41) 
(4.42) 



r 



One finds that i?i — > R + and i?2 — > -R- in the limit 
f — > 0, and -R 2 — > oo and R§ — > r as f — > oo. It 
is notable that Ri and i?2 are not independent but 
fulfills the constraint 



R 



H 1 U 2 



R 2 



(4.43) 



From R\ > R 2 > > -&l > -Rj, we conclude 
that if R 2 > Rj or R 2 < jR|, the ^ = constant 
curves are timclikc, while they are spacelike when 
R\ < R 2 < R\. In the region where R = constant 
curves are timelike (spacelike), a metric sphere is 
untrapped (trapped). 

We sketch in Figure [S] the region where R = 
constant curves are timelike (grey region), and 
spacclikc (white region). As r —> 0, we find that 
R± are null surfaces. 



C. Event horizons 

Let us proceed to demonstrate the structure of 
the future and past event horizons, utilizing several 
results obtained in the previous sections. 

Before embarking on this programme, let us re- 
capitulate basic features of event horizons. The fu- 
ture event horizon is defined by a future boundary 
of the causal past of future null infinity. The past 
event horizon is similarly defined by interchanging 
the role of future and past. These event horizons 
are by definition the achronal 3D null surfaces. In 
addition, the future (past) event horizon is gen- 
erated by null geodesic generators which have no 
future (past) endpoint in (M^g^) p2l |. 
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FIG. 5: Signature of R = constant surfaces for r = 1. 
R = constant curve is timelike (grey region) and space- 
like (white region), which are separated by trapping 
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horizons t 



If the spacetime is stationary, the black hole 
event horizon must be a Killing horizon [22j 
(see [49[ for the degenerate case). This theorem 
is considerably significant because it enables us to 
identify the locus of event horizon simply from the 
spacetime symmetry. It is a priori no interrela- 
tionship between the event horizon and the Killing 
horizon. Such a lucky consequence is exceptional 
to the stationary case. 

In the dynamical case, on the other hand, we 
have no specific guidelines for identifying the fu- 
ture event horizon but to evolve the spacetime into 
the infinite future. 4 Nevertheless, we can say, re- 



The Sultana-Dyer solution 13411 is an exceptional instance: 
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gardless of this adversity, that the black hole event 
horizon has to cover the trapped surfaces, provided 
the outside region of a black hole is sufficiently well- 
behaved [22||. Inspecting that the spacetime (|4.1[) 
appears to have a good behavior at least for f > 
and that the 9± — surfaces comprise null sur- 
faces (|4.32j) in the limit f — ► and t —> ±00, it may 
be reasonable to consider these null surfaces as the 
possible candidates of black and white hole hori- 
zons. Analyzing the near-horizon geometry and 
behaviors of null geodesies, we shall see below that 
this expectation is indeed true. 



the limit (|4.45p in Eq. (|4.ip . one obtains the near- 
horizon metric: 

d4 H =-T 2 ? 2 (l + t>)- 1/2 d? 

+ r-' 2 (l + tf) 1/2 (dr 2 + f 2 dft 2 ) , (4.46) 

which does not involve the parameter e. As a di- 
rect consequence of (|4.45|) . the near- horizon met- 
ric (|4.46[) is invariant under the flow 



t 



(4.47) 



1. Near horizon geometry 

From the behavior of trapped surfaces, we can 
deduce that the null surface f = is a plausible 
horizon candidate. We shall scrutinize the struc- 
ture of this surface in detail. 

We first look at the "throat" geometry (|2.8[) . 
Taking the fiducial time t = 1 for simplicity (the 
same conclusion is derived for any value of finite 
i), the proper distance s from the spacetime point 
(1, f , 9, <j)) to (1, r = 0,9, 4>) is given by 



s = lim / J g?f(l, f)df 



= lim 

f ^o 



r + In r 



(4.44) 



This implies that the point f = corresponds not 
to the regular origin of polar coordinates but to 
"spatial infinity," as in the extremal RN spacetime. 
In the extremal RN case, the future (past) event 
horizon is a null surface generated by f = and 
t = 00 (t = —00), with its infinite "throat" at f = 
with t being finite. Analogously, the event horizon 
of present spacetime, if it exists, should have r = 
with finite i as its "throat." So we deduce that 
the only candidate of event horizons in the present 
spacetime is f = and t = ±00. 

Therefore, the most convenient way to see the 
structure of these candidate horizons is to take the 
near-horizon limit, defined by 



t -> 



t 



er, 



0, 



(4.45) 



where e is a positive constant. Other conceivable 
limits fail to produce any sensible results. After 
the rescaling, the coordinate ranges of t and r are 
free from a restriction other than f > 0. Taking 



it is conformal to the Schwarzschild metric hence its causal 
picture is extracted in a simple fashion |5(J| . 



Namely, £ M is a Killing vector in the space- 
time P~36]) . 

Changing the coordinate f to the circumference 

radius R = (l + tf) 1 ^ 4 , the near-horizon met- 
ric (|4.46[) transforms into 



NH 



y 't 2 R 2 
16R 8 
(R 4 - l) 5 



R 5 



t(R A - 1) 



-dtdR 



-dR 2 



R 2 dQ 2 , 



(4.48) 



where 



f(R) :=t 2 (R 4 -l) 2 - R 4 

=t 2 (R 4 - R 4 + )(R 4 - Ri). 



(4.49) 



Here, R + and R- have been defined in Eq. (|4.32p . 
In the coordinates (|4T48|) . we have ^ = t(d/dt) fl . 
Apart from R = (which is indeed a curvature 
singularity) and the points at 9 — 0, n (which are 
north and south poles of 2-sphere), there appear 
additional coordinate singularities at t = and 
R = 1 in the metric ([Qg]) . 

Although the metric (|4.48[) is time-dependent, 
we can eliminate the time-dependence of the met- 
ric (|4.48p by changing to the time slice, 



ri ± := ln(±t) , for t ^ 0, 



(4.50) 



in terms of which the Killing vector is written as 
£ p = (<9/<977-j-) M and the near- horizon metric (|4.48l) 
is given by 



d5 2 

NH 



R? 
1Qt 2 R 1 

f(R) 



dm 



■dR 2 



4R 7 



-^dR 



(i? 4 - l)f(R) 



R 2 dtt 2 2 



(4.51) 



The sign of i]± has been chosen in such a way that 
n + (ry_) increases (decreases) as t increases. 

Performing a further coordinate transformation, 



AR 7 



(i? 4 - l)f(R) 



-^dR, 



(4.52) 
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the near-horizon metric (|4.51[) is brought into a fa- 
miliar form, 



dS 2 



NH 



f(R) 

R 2 



Wt 2 ^ 
f(R) 



di? 2 



(4.53) 



This metric describes a static black hole whose 
horizons occur where the lapse function vanishes 
f{R) = 0, i.e., where the Killing field ^ = 
{d/df±Y becomes null. The condition f(R) = 
gives two roots R = R± given by Eq. (|4.32[) . which 
coincide with the trapping horizons in the r — > 
limit taken. Thus, we conclude that the null sur- 
faces R = R± in the original spacetime are locally 
isometric to the Killing horizons in the static space- 
time (|4l53")) . 

Reminding the fact that the outside domain of 
the original spacetime (|4.ip is highly dynamical 
and hence is lack of non-spacclikc Killing field, it 
comes out a novel surprise for us that the near- 
horizon metric (|4.53[) permits the unexpected sym- 
metry (|4~47|) . Observe that the vector field (|3T4"T)) 
satisfies the Killing equation in the original space- 
time (|4.ip only at the horizon R = R±. This may 
be ascribed to the fact that the 11D solution is 
supcrsymmetric if all branes are at rest. The su- 
pcrsymmetry does not allow energy inflow, consis- 
tent with the property that the Killing horizon is 
totally geodesic. This may be clear by consider- 
ing the Raychaudhuri equation: T^ U ^S, V — > is 
indeed satisfied in the limit (|4.30j) or (|4.31|) . As 
far as the authors know, this is a first realization of 
asymptotic symmetry appearance at the black-hole 
horizon under the dynamical circumstance. 

Let us devote some space here to discuss the 
near-horizon static metric (|4.53|) in more detail. 
In this limit, the dilaton (|3.8[1 and the Maxwell 
fields (13.91) are reduced to 



i<5> = V6\nR, 



and 



K F [1 > = -4V2vrTir dT A di? , 
K p(S) = _ 4v /2^T J R 3 dT A di? . 



(4.54) 



(4.55) 



We can confirm Eqs. (|4"33")l . (|4"Jj4]) and (|4"35]) still 
satisfy the original field equations (|3.3|) . (|3.4[) and 
(|3.5[) . which justifies that the near-horizon limit 
(|4.45|) is well-defined. It is obvious that the near- 
horizon metric (|4.53[) describes a static black hole, 
whose asymptotic structure is neither flat nor AdS. 
Such an unusual asymptotic structure is, how- 
ever, that one commonly encounters in Einstcin- 
Maxwell-dilaton gravity (see e.g., [Hi], [III)- Albeit 
this peculiar asymptotics, it is easy to find that the 
causal structure is akin to that of the nonextremal 
RN-AdS solution (see Figure H]). The spatial in- 
finity R — > oo consists of a timclikc boundary X, 




t = oo 



FIG. 6: Conformal diagram of a maximally extended 
near-horizon metric (|4,53|) , The black hole has outer 
and jnner horizons, whose radii are respectively given 
by i?+ and i?_. Infinity consists of a timelike surface 
denoted by X. The white circles mark points at in- 
finity (i°, i or the "throat") and should not be re- 
garded as spacetime events. The filled circles B repre- 
sent bifurcation surfaces for the metric (|4,53[) at which 
£ M = (d/dT±) 11 vanishes, i.e., B's are the fixed points 
under this isometry. Thick red lines correspond to 
t — constant surfaces, and blue dotted lines denote 
f = constant surfaces. R = 1 is a coordinate singular- 
ity for the metric (|4.48p corresponding to t = 0. The 
shaded regions approximate our dynamical metric (|4.1|) 
in the neighborhood of horizons. 



a timelike singularity resides at the center R = 0, 
and the two distinct outer and inner horizons R± 
of black hole and white hole arise. 

Since i?+ is strictly larger than i?_ for finite 
t [see Eq. (|4.33[) ]. we find that the horizons are 
not degenerate. From the general formula k\ = 
=F(l/2)(V A1 ^ J/ )V' i ^ ly , one obtains the surface gravi- 
ties of these horizons, 



K± = ± 



f'(R±) 

8tR 5 , 



VI + 4r 2 
Vl + 4r 2 ± 1 



(4.56) 



4 It deserves to mention that the values of surface gravities 
are sensitive to the norm of the generator £ M of the Killing 
horizon. The spacetime H4.53H is not asymptotically flat, 
hence there exists no meaningful way to fix the values. (In 
the asymptotically fiat case, we usually require at infinity 



-1, which eliminates the ambiguity). However, 



this ambiguity causes no harm in our present discussion. 
The important point to note here is that they take non- 
vanishing positive values. We will revisit this issue when 
we discuss thermodynamics in the ensuing section. 
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where n + and k_ are the surface gravity of the 
outer and inner event horizons, respectively. The 
surface gravity of the Killing horizon is constant 
over the horizon, illustrating the equilibrium state. 
The nonvanishing surface gravity might be seem- 
ingly puzzling, in light of the situation that the 
11D solution was "maximally charged" because of 
the supersymmetry in the static limit. 
Introduce the null coordinates u± by 



= T± 
= In 



R Ar tdR 



f(R) 



±t\R 4 - irw 



R 2 + \^^\R 2 -R 2 _ l 1 ^- 
(4.57) 



where the integration follows from a direct calcu- 
lation by Eqs. (@3H|), and (|436) . The coor- 
dinates u± arc well-defined at R = Using u±, 
the metric (|4.53[) is translated into the single null 
form, 



dsl 



NH 



f(R) 

R 2 



du 2 ± - 2du±dR + R 2 dfll . (4.58) 



Killing horizon is the bifurcate Killing horizon. In 
their proof, it is assumed that the horizon is smooth 
(C°°-class), while the horizon in the present case 
is only finite times differentiable (C fe -class with k 
finite) (see Section TlV C 31 below) . 

To summarize this subsection, wc can expect 
that the original spacetime (|4.ip would have the 
future-event horizon at f — > with t — > oo, and 
the past-event horizon at f — > with t — > — oo. We 
have found that the point at f = with t being 
finite corresponds to "throat infinity" just as that 
of the extreme RN spacetime, at which future and 
past event horizons should intersect. In the neigh- 
, borhood of these horizon candidates, the space- 
time (|4.1[) is approximated by the near-horizon ge- 
ometry (|4.53p with Killing horizons, in which sev- 
eral portions of Killing horizons with radii R± ap- 
pear. What portion of Killing horizons in Figure |6] 
does it correspond to the horizon in our original 
spacetime? The answer is obvious: the "white hole 
portion" [grey-colored line segment encompassing 
blocks I, II and IIP in Figure [B] only satisfies the 
above criteria. 



Let us consider the plus-coordinate in Eq. (|4.58l) 
and discuss the outer white hole horizon (the 
boundary of I' and II). It is immediate to find 
that the null generator of the Killing horizon £ M 
is expressed in this coordinates as £ M = (d/du + ) fJ -. 
u + is the Killing parameter (£ M V fl u + = 1) of null 
geodesic generators. It then follows that the renor- 
malizcd tangent vector 



K< 



d 



(4.59) 



is affinely parametrized, i.e., fc M satisfies the zero- 
acceleration geodesic equation k"V v k^ = on the 
horizon. (Note that we are considering a white hole 
so that = This means that the 

affinc parameter A [fc M = (d/dX)^} is related to the 
Killing parameter u + as 



A 



(4.60) 



This manifests that the affine parameter A values 
from — oo to as u + ranges from — oo to oo, im- 
plying the bifurcation surface-a closed surface at 
which £ M vanishes. This can be verified by noticing 
^ = (d/df+) = n+\{d/d\y -> as A -> 0. The 
similar argument goes through to it_. Therefore, 
the nondegenerate Killing horizon (i.e., a Killing 
horizon with nonvanishing surface gravity) is in- 
complete cither into the past or future. This is a 
general consequence of a Killing horizon [2l| . 

However, the above discussion does not mean 
that the horizon in our original spacetime is the 
bifurcate Killing horizon. This is consistent with 
results in fl53j which asserts that the nondegenerate 



2. Null geodesies 

We marked out R + black hole horizon in the 
spacetime (|4.1|) . To conclude this more rigorously, 
we face up the problem of solving geodesic motions. 
Since the present spacetime is spherically symmet- 
ric, it suffices us to focus on radial null geodesies to 
argue causal structures. Although examinations of 
nonradial and/or timelike geodesic motions are im- 
portant issues in order to clarify the detailed phys- 
ical properties of the solution, we will not discuss 
these since behaviors of radial null geodesies are 
sufficient to determine the causal structure. 

The radial null geodesic equations are governed 

by 



l 



4H t ' ' 2f 2 H T H s 
t 2 {H s 



H 3 

—^f 2 =0, (4.61) 

4-H 

1 



\r 2 U 2 U% 



2H T 



-if 



4f 2 H T H s 



0, (4.62) 



-t 2 P + H T Hlr 2 = , (4.63) 



where the dot denotes a differentiation with respect 
to an affine parameter A. These equations are com- 
bined to give 



l ± r(H s + 3H T ) p 



2f 2 (ijf #5)1/2 



Cgrgl) 
2tH 7 



3U/2 



-2 



0. 



0, 



(4.64) 
(4.65) 
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where the plus (minis) sign refers to the outgoing 
(ingoing) geodesies. Unfortunately, the radial null 
geodesies do not appear to admit a first integral 
other than Eq. ()4.63|) . so it is not amendable to 
analytic study. Instead, we try to solve numeri- 
cally Eqs. (|4.64|) and (|4.65p subjected to the ini- 
tial constraint (|4.63j) . Making use of the degrees 
of freedom of the affine parameter A — > aX + b, we 
are able to choose A = at the starting point of 
the geodesies and set <(0) at any values we wish. 
Fixing the orientation of future-directed geodesies 

to be t > 0, and past-directed to be t < 0, we 
choose t(0) = 1 (—1) for future (past) directed 
radial null geodesies without loss of generality. 
Hence the residual freedoms that distinguish differ- 
ent geodesies are two, corresponding to the initial 
values [t(0) and f(0)] for each r > 0. 



Let us begin our consideration by the geodesies 
in the outside region f > 0. Taking the represen- 
tative spacetime events pi (I = 1,2,3) such that 

ti > t^j, < t 2 < £pjj an d t 3 < t^l where 
ti = t\ Pl (see Figure [3]), we have examined behav- 
iors of geodesies starting from t = tj. We call the 
geodesies emanating from the event pi as Class- 
I. Since fii > R + (R 2 < R-), Class-1 (Class- 
3) geodesies initially have a circumference radius 
larger (smaller) than R + (R-)- 

We depict several typical geodesic curves ema- 
nating from (i/,f(0) = 1) for r = 1 in Figure [71 
This is a representative figure for r < T cr j t . Quali- 
tative behavior of geodesies seems not so sensitive 
to the initial radial position f(0). The numerical 
results are summarized as follows: 



Class 1 Class 2 Class 3 

R R R 




FIG. 7: Radial null geodesies in the outside region (f > 0) emanating from t(0) = ti and f(0) = 1 for r = 1. 
The diagrams in the top, upper middle, lower middle, bottom rows correspond to future-directed ingoing null 
geodesies, future-directed outgoing null geodesies, past-directed ingoing null geodesies and past-directed outgoing 
null geodesies. The red and blue lines denote R+ and R-, respectively. 
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• Future- directed ingoing null geodesies: Class 1 
geodesies monotonically decrease the circumference 
radius and arrive at i?+ within a finite affine time. 
Class 2 geodesies first increase circumference radius 
since they are originally in the trapped region > 
0. But they always move across and finally 
reach R+ with decreasing area. The qualitative 
behavior of Class 3 are the same as that of Class 
2, except that Class 3 geodesies always cross R = 
R + twice, and if t-j is sufficiently close to t s , they 
may cross i?_ . All classes of geodesies have infinite 
redshift t(X) — > +oo when they finally arrive at 

• Future- directed outgoing null geodesies: Class 1 
geodesies necessarily go out to infinity R — > oo. 
Class 2 geodesies may extend out to infinity or ar- 
rive at the singularity t = t s if ti is small. Class 
3 geodesies inevitably plunge into the singularity 
t = t a within a finite affine time. 

• Past- directed ingoing null geodesies: Class 1 



geodesies originate from R > R + and their 
radii monotonically decrease toward Class 
2 geodesies have qualitatively the same behavior. 
Class 3 geodesies start from R < 7?_ with increas- 
ing area, then cross R- (with finite t), attain the 
maximum radius and get back to i?_ again with 
undergoing infinite blueshift. 

• Past- directed outgoing null geodesies: Class 1 
geodesies may initially increase the area, but all 
geodesies unavoidably terminate into the singular- 
ity t = t s within a finite affine time. 

From these results, we conclude that the null 
surfaces R = R± locate within a finite affine time 
from outside spacetime events. Behaviors of future- 
directed outgoing null rays of Class 2 geodesies im- 
ply that there exists a critical null curve t = i*(f) 
such that outgoing rays emanating from t > t * can 
get to infinity, whereas outgoing rays emanating 
from i < fall into the singularity. 



Let us discuss next the geodesies inside the hori- 
zon. We call Class 1 as t > i\ > t^Hi Class 2 
as > t 2 > and Class 3 as > t 3 (see 
Figure U]) . Figure H] plots the geodesic curves ema- 
nating from the spacetime event (tj, f (0) = —1/10) 

with t(0) = ±1. Geodesies starting from f (0) < r 
show the same behavior as Class 1. The result is: 



• Past-directed outgoing null geodesies: All 
geodesies eventually approach R = R + with infi- 
nite blueshift t(A) — > — oo. 

From these results, the null surface R = R_ has 
an ingoing null structure, analogous to the white 
hole horizon or the black hole inner horizon. 



• Future-directed ingoing null geodesies: Class 1 
geodesies and Class 2 geodesies starting from not so 
small t2 eventually fall into the singularity t = t s . 
Class 2 geodesies with t 2 close to t^(f) fall into 
the singularity r = — 1 within a finite affine time. 
Class 3 geodesies initially increase area, but they 
eventually plunge into the singularity at f = — 1. 

• Future- directed outgoing null geodesies: Class 1 
geodesies increase area and approach R = R- with 
infinite redshift t(X) — > +oo. Class 2 and Class 3 
geodesies decrease area and arrive at R = R_ with 
infinite redshift i(A) — > +oo. 

• Past- directed ingoing null geodesies: Class 1 
geodesies monotonically reduce area and fall into 
the singularity t = t s . Class 2 and Class 3 geodesies 
initially grow the area and then turn to decrease, 
and finally fall into the singularity f = — 1. 



3. Asymptotic solutions of geodesies 

We have numerically established that the ra- 
dial null geodesies are incomplete at the null sur- 
faces R = R±. We shall look into the asymptotic 
geodesic solutions and discuss further the horizon 
structure. 

If a null geodesic is known as t = i(r), the affine 
parameter A is obtained by a simple quadrature [38| 



A = / df exp 



± / U[t(r'),r']dr' 



(4.66) 



(4.67) 



where we have used the shorthand notation 
2tH t ' 

In the event horizon limit (|4.30p . we have 

U -> — . (4.68) 



f(l + Vl + 4r 2 ) 
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Class 1 Class 2 Class 3 




R R R 




R R R 




FIG. 8: Radial null geodesies in the inside region (f < 0) emanating from t(0) — ti and f(0) = —1/10 for r = 1. 
The diagrams in the top, upper middle, lower middle, bottom rows correspond to future-directed ingoing null 
geodesies, future-directed outgoing null geodesies, past-directed ingoing null geodesies and past-directed outgoing 
null geodesies. The red and blue lines denote R+ and R-, respectively. 



Substituting this into Eq. (|4.66l) and solving with 
respect to the radial coordinate, we obtain the 
asymptotic solution of the future-directed null 
geodesic around the horizon R = R + as 

r = c{ + \\ - , t = 4 +) (A - A+r 1 /^ , 

(4.69) 

where A+ corresponds to the arrival time for 
the geodesies at the horizon, and and 
are positive constants satisfying c^c^ = (1 + 
Vl + 4t 2 )/(2t 2 ). k + has been given in Eq. (|4.56j) . 
We can find from Eq. (|4.69j) that the radial 
geodesies indeed reach the horizon within a finite 
affine time [HI . Equation (|4.69|) implies that f and 
t are not smooth functions of A [note that 
never takes an integer]. 
Similarly, we obtain 

r = C [-\X - X.) 1 ^- , t = c[-\\ - A.)" 1 /- , 

(4.70) 

for an outgoing null geodesic near the horizon R = 



R-. Constants c\ and c 2 satisfy c\ c 2 = 
(1 - \/IT4^)/(2t 2 ). 



D. Carter-Penrose diagram 

We are now in a position to discuss global causal 
structures of spacetime, by assembling considera- 
tions hitherto obtained. The optimal way to ap- 
preciate the large scale causal structure is to draw 
the Carter-Penrose conformal diagram, which en- 
ables us to visually capture the global light-cone 
fabric. We first notice the followings: 



(i) The only candidate of future and past event 
horizons are f = and t — > ±oo, which are 
joined at the "throat" at r = and t being 
finite. 

(ii) The near-horizon geometry of the event hori- 
zons is locally isometric to that of the static 
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black hole (|4.53[) . The white- hole portion cor- 
responds to the horizon in the original space- 
time (|2.1j) with the bifurcation surface re- 
placed by a smooth surface. 

(iii) There are the curvature singularities at t = 
t s (r) = — 1/f and r = — 1. These singular- 
ities are timelike (section lIV Ap . The time- 
dependent singularity t s (r) present in the 
f < domain exists for t > 1, whereas t s (f) 
lying in the f > region and the other sin- 
gularity f = — 1 exist for eternity. 

These observations prompt us to imagine the po- 
sitional relation between singularities and the hori- 
zons. Figure H] describes the conformal diagram 
of our dynamical black hole. From properties (i) 
and (ii), we can depict the identical horizon struc- 
ture as in Figure [S] Since the t = constant(< oo) 
lines are everywhere spacelike, each slice originates 
from the "throat" f = 0. For negative values of t, 
t = constant surfaces must intersect the singularity 
t s (f) at finite r (see Figure [3]). Considering prop- 
erty (iii) that the singularity outside the horizon 
is only t = t s (f) < 0, the right side dotted por- 



tion of grey line can be drawn in Figure [3] Outside 
the horizon r > 0, one can depict the contours of 
t = constant and f = constant family of surfaces, 
both of them to be orthogonal (Figure [5]). These 
aspects are all consistent with our numerical survey 
of geodesies. We thus conclude that the spacctime 
metric (|4.ip indeed describes a black hole in the 
FLRW universe [aside from the undesirable time- 
like naked singularity t s (f)]. Although the null sur- 
face R- is a one way membrane of "region of no 
entrance," it does not deserve to be a white hole 
horizon in a mathematical sense since the space- 
time does not possess the past null infinity. 

Inside the event horizon R < R + , the timelike 
singularities are vertically joined at t = — 1. The 
past boundary R + can be matched to the black hole 
horizon. We can find as sketched in Figure |H] that 
these patches are infinitely arrayed vertically. It 
should be emphasized, however, that this is only a 
possible extension, since the horizon is not analytic 
in the present case. One may glue the near-horizon 
geometry (|4.46[) to the spacctime (|4.ip across the 
horizon. 




(a) The contour of t = constant. (b) The contour of f = constant. 



FIG. 9: Conformal diagram of the spacetime (|4.ip . We draw (a) the contour (spacelike) curves of t = constant 
by red lines and (b) the contour (timelike) curves of f = constant by blue lines. The singularity of R — consists 
of three parts: a black- hole singularity at f = — 1 and at f = r s (t) = — 1/t inside the horizon, and a "big bang 
singularity" t = t s (r) = —1/f outside the horizon. 



Next, we wish to fill in the trapping horizons into (i) There arc trapping horizons t — at which 

this diagram (we only consider the r = lease). We § T = fj. Outside the horizon (f > 0), the 

should remind the following remarks: , n , • c, i • r(+) • . 

6 whole portion of trapping horizon is past- 



22 



outer for r < T cr i t hence always spacelike, 
analogous to FLRW universe filled by a stiff 
matter. While the trapping horizon is 
always timelike. Inside the trapping horizon 
(f < 0), £™ coincide at r = fo- i^H anc ^ a 
part of near r ~ ro are timelike. Other 

portion of f = changes signature near 
f = into spacelike. 

(ii) Trapping horizons occur where R = constant 
surfaces becomes null. The contour curve of 
circumference radius is spacelike for R2 < 
R < R\ (see Figure [5]) . As approaching the 
event horizon, the trapping horizons tend 
to have constant radii R±. 

(iii) For f < 0, the circumference radius R = 
[(I + tr)(l + f) 3 ] 1 / 4 becomes infinitely large 
as f -> —00 with f staying constant. We 
can show following the same argument in 
Eq. (|4.44|) that this is a "past timelike in- 
finity." 5 

Outside the horizon f > 0, the R = constant (> 
R+) surfaces are the same as the FLRW cosmol- 
ogy: there exist a spacclikc past trapping horizon 
£™(f), above which R = constant surfaces are 
timelike and below which R = constant surfaces are 
spacelike (see Figure. [10]). For -Jt = 1 < R < R + . 
R = constant curves are everywhere spacelike and 
lie in the future of a critical null curve t = 
For iJ_ < R < 1 = y/r, R = constant curves cross 
the future trapping horizon and change sig- 

nature. For R < R = constant curves are 
always timelike. 

Inside the horizon (f < 0), R = constant curves 
are the same as outside for R < i?_ and ^fr = 1 < 
R < Whereas, R = constant curves for R > 
R+ quite differ from those in the outside. They 
cross the trapping horizons twice. 

V. BLACK HOLE THERMODYNAMICS 

Black hole thermodynamics has been established 
as rigorous mathematical laws of black holes with 
Killing horizons [231 - 4301 . Since the three laws of 
black hole thermodynamics interrelate the clas- 
sical gravity, quantum mechanics and statistical 



5 This is slightly different from the extremal RN geome- 
try, for which it takes an infinite affine time to reach the 
corresponding point from inside the black hole, but the 
point locates at the finite circumference radius. The rea- 
son why the point in the present spacctime has an infinite 
circumference radius inside the horizon might be due to 
the cosmic expansion. 




FIG. 10: A conformal diagram of the black hole in 
the expanding universe. The curves of R = constant 
are shown by green lines. Null curve _R+ is the future 
event horizon, and T + is the future null infinity. The 
trapping horizons t = t^(r) and t = t^H^) are a ^ so 
shown by red and blue lines. The diagram can be ex- 
tended beyond the upper i?_ null curve (possibly with 
the contracting patch r < that turns the above figure 
upside-down) in a continuous but a nonanalytic man- 
ner. 

mechanics, they are likely to have a key role to- 
ward quantum laws of gravity. Lastly, we discuss 
the thermodynamic properties of the present time- 
dependent black hole. 

Since the present spacctime (|4.1[) possesses a 
Killing horizon, the thermodynamic laws continue 
to hold without change. It turned out that the 
Killing horizon has nonvanishing surface gravi- 
ties (|4.56p . However, their values are dependent 
on the normalization of the Killing field generators 
of the horizon. Here, we intend to obtain the tem- 
perature associated to the time translation in the 
FLRW universe. To this end, we resort to the laws 
of trapping horizons. This is a hotly-discussed is- 
sue in recent years [3 EE 113 ■ Laws of trapping 
horizons are local extensions of black hole thermo- 
dynamics. 

When discussing the dynamical aspects of black 
holes, a major obstacle for extending the black hole 
thermodynamics to a non-stationary setting is how 
to define a surface gravity. In the case of a spheri- 
cally symmetric spacctime, a natural time direction 
is specified by the Kodama flow [3, HE [H| . 

Write the spherically symmetric metric as, 

ds 2 = g AB (x)dx A dx B + R 2 (x)dnl , (5.1) 

where qab (x) (A, B = 1, 2) is the two-dimensional 
Lorentz manifold (M 2 , gAB) perpendicular to the 
metric sphere. The coordinate, x A , on M 2 corre- 
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sponds to t and r. Using this coordinate patch, the 
Kodama vector is defined by [54| 



k a = _ € ab DbRj 



(5.2) 



where sab and Da are the volume element and 
the covariant derivative of qab- We may view 
K A as a spacctime vector by = K a (8aY '■ ^ 
follows immediately from the orthogonal property 
K^^R = that is divergence-free, V„K» = 
R~ 2 D A (R~ 2 K A ) = 0. Another key property 
comes from the relation K^K^ = -(V fJl R){V fl R), 
so that is timelike (spacelike) in the untrapped 
(trapped) region, i.e., defines a preferred time- 
like direction in the untrapped region, irrespective 
of the non-stationarity of spacctime. Specifically, 
the Kodama vector becomes null at the trapping 
horizon, just as in the same way the Killing vector 
becomes null at the Killing horizon. 

It is enlightening here to look into the rela- 
tion between K^ L and the Misner-Sharp energy, 
which is also characteristic to spherically sym- 
metric spacetimes [43|. Inspecting V M if„ = 
DaKb(V p,x A )(Vv% B ), a simple calculation shows 
that Q nv^J ^ K v = holds in any spherical space- 
timcs. Hence we can define a divergence-free vec- 
tor field k 2 J M = —Q^ W K V , representing an energy 
current due to the Einstein equations. Integration 
of over the volume V with exterior boundary S 
yields the Misner-Sharp energy m = — j v J^dS^. 
To summarize, the Misncr- Sharp energy is a charge 
associated with the locally conserved current. 

As seen above, the Kodama vector in the spher- 
ical spacctime plays a role similar to the Killing 
field in stationary spacetime. One can speculate 
that laws of trapping horizons are related to an 
observer along the Kodama flow. 



irpjj are given by 



,-.\3 



(+) _ (l + f) 6 (l+4f) 



8i?f 



4f 2 + 8r+l Bs 



(5.4) 



_ ( _) _ (l + f) 6 (l+4f) 



4f 2 + 8f + 1 - 



m-i 



(5.5) 



Taking the event horizon limit f — > [see 
Eqs. (|4.30|) and (|4.31[) ] in the above equations, we 
obtain the black hole temperature 



r, 



(±> 



TH 



BH 



2nQ 



VI + 4r 2 
16ttt 2 R\Q 



(5.6) 



Comparing this with Eq. (|4.56p . these are equiva- 
lent to surface gravities associated with renormal- 
ized generator of the horizon, 



( d 



\df ± 



1 / d 



4tR 3 , \df ± 



(5.7) 



which coincides with the Kodama vector evaluated 
on the horizon for the near- horizon metric (|4.53l) . 
For the future horizon the temperature 



BH 



takes the maximum value 



(+) 



BH(max) 



2irQ 



0.00213856 Q- 



(5.8) 



at t = a/3/2 for fixed charge. The temperature 
vanishes in both limits of r — > oo (degenerate hori- 
zon) and t — > (no horizon). The former recovers 
the result for the extremal RN black hole. The 
temperature at the past horizon on the other 
hand, has no maximum value. It monotonically in- 
creases to infinity as r — > 0. is always higher 



than T, 



(+) 



BH 



A. Temperature: Oth law 

A naive definition of the surface gravity for the 
trapping horizon is to replace the Killing field by 
the Kodama vector in the definition of surface grav- 
ity of a Killing horizon. This prescription does not 
work, since the trapping horizon is not the null sur- 
face generated by the Kodama vector. A proposed 
definition of surface gravity for the trapping hori- 
zon [44|, |4{| is given by the "equilibrium part" 



(5.3) 



where the equality is evaluated at the trapping 
horizon. After some amount of algebra, one finds 
that the surface gravities of trapping horizons t = 



B. Energy balance: 1st law 

It is a widely accepted criterion that a well- 
defined energy should satisfy an energy balance 
law. The Misner-Sharp energy indeed fulfills this, 
as in Eq. (|4.17[) . We can rewrite each term in this 
equation into more recognizable form. Defining 2D 
quantities, 



Ya = t ab d b r + pd a r, 

(5.9) 



and making use of the Einstein equations, one ar- 
rives at the unified first law [45| , 



D A m = Aip A + P c ffD A V . 



(5.10) 
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where A = 4ttR 2 and V = ^j-R 3 denote the area 
and volume of the metric sphere. This equation 
illustrates that the mass variation is supplied by an 
injection of energy current and the external work 
term. The expression of ipA is comparatively messy, 
but it is straightforward to obtain. 

Projecting Eq. (|5.10j) along the generator , C,^ = 
C j4 (9 j 4) m , of the trapping horizon and noticing the 
fact that C /i V A1 (m/i?) = 0, one obtains 

At A ip A = ^( A D A A. (5.11) 

This is an energy balance law of a trapping horizon 
(see [H, Hi| for detailed derivation). Here, along 
the trapping horizon with 6+ = 0, £ A is obtained as 
Q A = ±e AB Db0+ (where the sign should be appro- 
priately chosen in such a way that it is outgoing in 
the spacelike case or future-directed in the timelikc 
case). 

Using Einstein's equations, we find that the sur- 
face gravity is expressed in terms of the Misner- 
Sharp energy and the pressure as 

K ™ = ^ - 4 ^ p cff , (5.12) 

where equality is understood at the trapping hori- 
zon. This is the Newtonian analogue definition of 
acceleration. 



C. Entropy: 2nd law 

It follows from the first law of a trapping horizon 
that we can identify the entropy by one quarter of 
the area of the trapping horizon, i.e., it accords 
with the "Bckcnstcin-Hawking formula," 

a(+) _A 1 _ irRl ( _) _ A 2 _ txR\ 

— AG ~ G H _ 4G _ G 

(5.13) 

Taking the event horizon limit (14.301) or (|4.31[) . we 
recover the well known result [2f| |27[ 

In the limit t — > oo with fixed charge, the above 
entropy reduces to that of the extremal RN black 
hole, S^h = ttQ 2 /G. 

VI. CONCLUDING REMARKS 

In this paper, we have made a thorough discus- 
sion about the causal structure and physical prop- 
erties of the spacetime derived from intersecting 
M-branes. We have found that the solution in- 
deed describes a black hole embedded in the FLRW 



cosmology filled with fluid obeying the stiff equa- 
tion of state. The global causal structure is dis- 
played in Figures and [101 Since the solution is 
approximated by the extreme RN solution near the 
"throat" and the flat FLRW universe with P = p 
at infinity, one might first envisage that the causal 
structure is obtainable by patching these two lim- 
iting spacetimes. That is to say, according to our 
first intuition, one might have expected that space- 
time should possess a spacelike big-bang singular- 
ity at t = 0, there should exist a degenerate event 
horizon, and the timelike singularity should appear 
only inside the hole. However, our careful analysis 
revealed that the global causal structure is com- 
pletely different from the above rough estimate. 

Our solution satisfies the dominant energy con- 
dition, so that the energy densities are always posi- 
tive, absolute value of principal pressures do not ex- 
ceed the energy density for respective fluids and the 
energy flux current is always causal. This desirable 
property is not seen in the solutions found in the lit- 
erature. Hence the results presented in this paper 
open up new avenues for further research on black 
holes surrounded by usual matters in the expand- 
ing universe from higher-dimensional point of view. 
Our solution, however, may not have a direct astro- 
physical relevance because of nonzero electromag- 
netic charge. The charge is probably also responsi- 
ble for the timelikc singularity t — t s (r) outside the 
horizon. The timelike singularity does not develop 
as the big-bang singularity in the usual FLRW cos- 
mology with fluid P = wp (—1 < w < 1). Unfortu- 
nately, the construction of black hole solution with- 
out charge may be beyond the intersecting brane 
picture. 

In the process for obtaining the global structure, 
we gave a coherent description concerning the trap- 
ping horizons. The main idea on which our dis- 
cussion is based is that the trapping horizons re- 
flect the physical situations of marginal surfaces 
on which the either of expansions of light ray van- 
ishes. This local character enables us to relate it 
to the curvature singularity and the Misner-Sharp 
energy. A more important belief to which we resort 
is that the trapping horizon with negative outgoing 
expansion does not occur outside the horizon. The 
present spacetime indeed has this property (except 
in the neighborhood of singularity). We confirm 
the infinite redshift (blueshift) surface as a black 
hole horizon (white hole horizon in a quoted sense) 
combining the analysis of near-horizon geometry 
and the behaviors of null geodesies. 

It was somewhat surprising that the solution ad- 
mits a nondegenerate Killing horizon. The Killing 
horizon is usually associated with symmetry of 
time-translation and angular-rotation. The black 
hole remains the same size and fails to grow, al- 
though the the black hole is surrounded with fluid. 
This characteristic property may be ascribed to the 



25 



fact that the 11D solution was supersymmetric in 
the static limit. Although the dynamically brane 
intersecting solution breaks supersymmetry, it still 
maintains a part of the BPS characters. The same 
takes place in the Kastor-Traschen black hole. 

In this paper, we have taken a particular notice 
on the solution, whose 11D "oxidized" solution has 
four kinds of harmonics of spherical symmetry. A 
more general non-spherical spacetime is of course 
more complicated. Still, the profound understand- 
ing of spherically symmetric case will of substantial 
aid in exposing more complex structures of dynam- 
ical black holes. 
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Appendix A: Intersecting brane and black 
holes 

In this Appendix, we consider an intersecting 
brane system in the 11D supcrgravity theory, which 
is expected to be an effective field theory of M- 
thcory. We discuss how to obtain the 4D effec- 
tive action and produce solutions in the 4D space- 
time. We intend to consider M-branes, for which 
the Chern-Simons term " F A F A A" has no con- 
tribution. Hence, it suffices in our setup to concen- 
trate on the following 11D effective action, 



S = 



2k 11 



d 11 A v ^T 



1 



2(p A + 2)\ 



{Fp A +2) 2 



(Al) 



Here, A denotes the type of branes with which the Abelian (pa + 2)-form field F PA+ 2 is coupled, and 
Pa (=2 or 5) is the dimensions of branes. 

Once the 11D brane configuration is given, the 4D solution is derivable via the standard toroidal com- 
pactification. We shall analyze intersecting brane systems involving four-charges, which is needed to find 
a 4D maximally charged (supersymmetric or nonsupersymmetric) black hole with regular event horizon. 
We can construct two kinds of such a configuration: M2-M2-M5-M5 and M2-M5-W-KK system. As a 
concrete example, we focus our attention to the M2-M2-M5-M5 intersecting brane system. Compatibility 
with 11D supcrgravity equations of motion determines respective brane codimensions, which is given by 
Table U 
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10 


M5 
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o 












M5 
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o 








M2 


o 








o 




o 










M2 


o 










o 




o 









TABLE I: M2-M2-M5-M5 brane system, 
branes. 



The circles describe which dimensions are filled by the corresponding 



There appear four charges (Q2, Q2', Q5, Q5') associated with the corresponding four branes. In the 
static spacetime, we have the following intersecting brane solution: 

+H^H^dy 2 + H^H^dy 2 + H^H^dy 2 + H^H^dy 2 + (dr 2 + r 2 d^)] , (A2) 

I 

where Ha are harmonics on the three Euclidean the metric in the square bracket correspond to the 
space (ds 2 = dr 2 + r 2 d£l 2 ). One can see immedi- dimensions to which the A- brane belong, 
ately that the directions involving inverse Ha for 



2G 



If one toroidally compactifies common 7D world- 
volume of branes, a 4D solution is obtained. 
Rewriting the 11D metric as 

7 7 



where 



b\ = b\ = b\ 



2 \ 1/3 



H2H21 
H5H51 



1/3 



(A4) 



1 "5 



b 2 - 

6 _ 1 flllfo 



H 2 ,HlV /3 



, hi, 



2 \ i/3 



,(A5) 



and compactifying y^-coordinatcs (i = 1, ••■ ,7), 
the 4D solution in the Einstein-frame is given by 



ds\ = -Edt 2 + - (dr 2 + r 2 dft 2 



2 1 



where 



-1/2 



(A6) 



(A7) 



If we assume that harmonics Ha's arc spherically 
symmetric, i.e., 



Ho = 1 



Hk = 1 



H 2 , = 1 + 9*. , Hw = l + ^. (A8) 
r r 

we find a static extreme black hole solution in 
four dimensions. Here, Qa"° represent the branc 



charges. If all charges vanish, both of the 11D and 
4D solutions are trivial. The extension of the har- 
monic functions Ha as discussed in Appendix |B] 
gives a multi-black hole system. 

Generalizing this solution to the time-dependent 
one, we find that only one brane among four can be 
time-dependent under the metric ansatz assumed 
in [20I ] . The intersecting brane metric is still given 
by Eq. (|A2j) . The field equations require that the 
time-dependence is linear, i.e., the metric functions 
in the spherically symmetric case are given by 



Hn 



H s 



1 + 9l 

t r 
r 



1 



H s 
H s » 



^5 

1 

r 



(A9) 



where to is a constant with dimension of time. 
T and S,S',S" denote one time-dependent brane 
and three static branes, respectively. Any one 
of M2, M2', M5, and M5' branes can have time- 
dependence. This gives a black hole in the expand- 
ing universe in four dimensions, which we discuss 
in this paper. It is also extended to a multi-black 
hole system (see Appendix [B]) . If all brane charges 
are set to zero, 11D solution is the Kasner solution 
describing a homogeneous but anisotropic vacuum 
universe, whereas the 4D solution reduces to the 
flat FLRW cosmology. 

Since we know the 11D action (|Aip . assuming 
the brane configuration shown in Table U and com- 
pactifying the spatial directions as (|A3|) , we can de- 
rive the effective 4D action, which gives the present 
time-dependent solutions, as follows. 

The scales of extra dimensions bi (i = I, • • • , 7) 
behave as scalar fields in 4D spacetime, i.e., the 
effective action of gravity sector is 



J 



Si{g,h) 



d 4 x 



V^lnM +2^(Vln6i 



(A10) 



i=i 



i=l 



Although we compactify seven dimensions, we have only four branes. Hence the degrees of freedom are 
maximally four. This can be affirmed by writing down the kinetic term of the scalar fields in terms of 
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harmonic functions Ha as 

\ fv^ln^ +2^(Vln6 l ) 2 

= [3{(VlniJ 2 ) 2 + (VlniJ 2 <) 2 + (V\nH 5 f + (V\nH 5 ,) 2 } -2(Vlnff 2 • VlniJ 2 ' + V In H 5 ■ V In H 5 >) 
- 2(Vlnff 2 ■ VlnJ? 5 + Vlnif 2 ■ VlniJg, + Vm# 2 , ■ Vln# 5 + Vlnff 2 ' ■ Vln# 5 ,)] 

= ^ [(Vln(H 2 /H 2 ,)) 2 + (Vln(H 5 /H 5 ,)) 2 + (V \n(H 2 / H 5 )) 2 
+ (W\n(H 2 /H 5 ,)) 2 + (W\n(H 2 ,/H 5 )) 2 + (V \n(H 2 , / H 5 ,)) 2 ] 

4e( v m 2 . (ah) 

A<B 



where 



$ survives, which is normalized from (|A11[) as 



K0 AB 



2V2 



In 



Ha 
H b 



(A12) 



K(pTS = 



V3 



111 



2^2 \H S 



H, 



denotes the "scalar field mixing" term. 

Supposed that all charges are equal (Q 2 = Q 2 i = 
Qh — Qb' = Q) as in the main text, it follows that 

4>SS',4>S'S",<t>SS" are trivial, and 4>ts,<I>ts' ,4>ts" 



This is identical to Eq. f|3T8|) . 

Next we reduce the 11D form- field sector as fol- 
lows: the M2 and M5 branes couple to four-form 



are the same. As a result only a single scalar field and its dual seven-form field, respectively. 
I 



Hence the effective action of the form fields is reduced to four dimensions as 

i^E5(^(w=(n».-') 2 n^»< (n».-) 3 



gWgV a 



r p,v r pa °4 °6 + e jiv r pa °5 °7 + jxv r pa °1 °2 °3 °4 °5 + /if r pcr °1 °2 °3 °6 °7 



= v-g X 



3/2 



(^) 2 + 3 



#7 



-1/2 



(A13) 



where we set 



2tt 



k 2 ( P a + 2)! 



PA+2J 



This ansatz is consistent with our result p. lip for 
the Maxwell fields, because the electric potential of 
four- form field in 11D is given by = I/Ha + 



(t), where aS A >(t) is an arbitrary function of t, 



which comes from a gauge freedom. As a result, 
we obtain the effective 4D action for the form-field 
sector as 



+3e 



±V6k*/3 



which is the same as Eq. 



In the static case with equal charges, the 4D so- 
lution (| A6|) with (|A8|) corresponds to an extreme 
(A14) RN black hole (which is indeed a solution in the 
Einstein-Maxwell system). While for the time- 
dependent case with equal charges (|A9[) . it de- 
scribes a black hole in the FLRW universe which we 
have established in the body of the present paper. 



Appendix B: Multi black holes in the 
time-dependent universe 



Writing Ht = t/to+Hx, Hs = 1+3$ and so on, 
the 11D supergravity equations of motion require 
the functions Ha {A = T, 5, S', S") to be arbitrary 
(A15) harmonics on flat three space. Hence, just by re- 
placing the monopole term Qa/t by multicenter 
harmonics, wc obtain a collection of black holes in 
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a dynamical background. To be specific, we have 

(Bl) 



N n (A) 

^ = Er 



(A), 



where the constants and Q\ A ' (> 0) corre- 
spond to the loci and the charges of i-th black hole 
associated with A-branes, respectively. The lin- 
ear term (aj • r%) has been dropped by the asymp- 
totic boundary conditions at infinity. In the case of 
three equal harmonics Hs = Hs' = Hs" , this met- 
ric solves the field equations of Einstein-Maxwell- 
dilaton system (|3.2p if the dilaton and U(l)-gauge 
potentials are given by (|3.8|) and p. lip . 

Near each mass point (with t being finite) r\ , 
there exists an infinite throat as in the single mass 
case discussed in the body of text. Far from the 
throat, on the other hand, the solution tends to an 
FLRW universe filled by a stiff matter. 

As in the case of the Kastor-Traschen solu- 
tion [37], [38[ , this spacetime is expected to describe 
a collision of black holes provided the background 
universe is contracting (to < 0). Since each black 
hole is ignorant of others, i.e., the gravitational and 
electromagnetic forces between black holes are bal- 
anced, the collision occurs by a brute-force method 
responsible for the background contracting uni- 
verse. The difference from the Kastor-Traschen 
case lies in that the background universe obeys the 
power-law contraction a oc t 1 / 3 , so our discussion 
parallels [l|| in which colliding D3 branes were dis- 
cussed in detail. 

Let us start with the negative time t < and 
run time forwards. Since the t = surface is again 
nonsingular, the universe continues to shrink for 
positive values of t until the singularity Ht = is 
reached. Specifically, the metric continues to ex- 
ist inside the domain, D t , bounded by the level 
set Ht = t/(— to). It then follows that at small 
positive t, the domain D t is a large connected vol- 
ume containing all black holes. As time passes, 
the domain D t continues to contract and tends to 
spilt into disconnected pieces containing each mass 

(A) 

point r\ . This means that black holes scatter off 
rather than coalesce, and the universe is bounded 
by curvature singularity at Ht = 0. 

The multiple black-hole solution (|B1[) is to be 
compared with the Kastor-Traschen solution, 



ds 2 



-H~ 2 dt 2 + H 2 dr 2 



with kF = d(iJ _1 ) A dt and 



H 



t N 



Qi 



to 



r - n 



(B2) 



(B3) 



This is an exact solution of Einstein-Maxwell-A(= 
3/fg) system. The distinction between our space- 
time and the Kastor-Traschen solution is essen- 
tially only the power of the lapse function, where 



divergence of lapse corresponds to the curvature 
singularity for each solution. The exponent is 
closely associated with the number of branes and 
more general class of solutions is available. Further 
detailed analysis will be reported elsewhere (55j . 



Appendix C: 5D time-dependent black holes 

As proved in [2(j, a 5D time-dependent "black 
hole solution" is also obtained from the M2-M2- 
M2 and M2-M5-W intersecting brane systems. Let 
us discuss the former case. As in the static coun- 
terparts, we need only three nontrivial charges to 
obtain a black hole solution. The 5D metric in the 
Einstein frame reads 



ds 2 



-S 2 dt 2 + (dr 2 + r 2 dn 2 ) , (CI) 
with 

Z = (H T H s H s ,y 
Here we have introduced 



-1/3 



(C2) 



ffj 1 - 



t , Q q 



to ' 



H s = 1 



Ho, = 1 



(C3) 



to denote the harmonics in the flat 4D space. Qt 
and Qs, Qs' are charges of one time-dependent and 
two remaining static M2 branes, respectively. The 
lapse function S takes a relatively simple form com- 
pared to the 4D metric (|2.2|) . 

Assuming t/to > and transforming to the new 
time coordinate i given by 



— = ( — V with tn = — - 
to Vo 2 



(C4) 



the solution (|C1J) is cast into the form, 



ds 2 = -E 2 dP + ^ (dr 2 + r 2 dfi|) , (C5) 



where 



(HtHsHs' 
(i/tof 4 , 



-1/3 



with 



Hn 



1 



Qi 



(C6) 
(C7) 

(C8) 



The expansion law with (|C7|I is again identical 
to that of the 5D universe with a stiff matter. The 
limit of r — > with keeping t finite gives the same 
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"throat" geometry of the 5D extreme RN black 
hole. According to the detailed argument laid out 
in the main text, we may regard this solution as a 
black hole in the expanding universe. 

For the case in which Qt = Qs = Qs' == Q, 
the 5D metric is an exact solution of Einstcin- 
Maxwell-dilaton system whose action is given by 



16tt 



if the dilaton is given by 



(C9) 



where 



(C14) 



R + and i?_ correspond to the future event hori- 
zon and the past event horizon, respectively. The 
Carter-Penrose diagram is quite similar to that in 
4D, although there exist minor differences. 
The surface gravities are found to be 



,(±) 

"BH 



\/T + 16t 2 
12r 2 Q 1 /2 J R7 



(C15) 



^ 1 , ( H T 



(CIO) 



and the electromagnetic fields take the form, 



±1 T 



with coupling constants 
4 



At = 



V3' 



A<? — A.s' — — 



V3' 



(Cll) 



(C12) 



After short calculations, we find the following 
results. The horizon radii are given by 



1? 



Q 



3/2 



+ 16t 2 ± 1) , 



(C13) 



The temperature of the future event horizon 

Pth = ^bh/X^ 71 ")] vanishes in the both limits of 
r — >• and of r — > oo, just as the 4D black hole. 
The maximum temperature is given by 



T, 



(+) 



4 Vl3 V 7 ! + VTV62 + UVU 



BH(max) 



tt(6+ V62~ 
0.0395465 Q 



^\ 7/3 



-1/2 



at 
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0.364381 . 
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